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O Abstract 

Li-. We study operators of Kramers-Fokker-Planck type in the semiclassical 

^^ limit, assuming that the exponent of the associated Maxwellian is a Morse 

'""' function with a finite number no of local minima. Under suitable additional 

^ assumptions, we show that the first no eigenvalues are real and exponentially 

small, and establish the complete semiclassical asymptotics for these eigen- 



%^ values. 



Resume 

Nous etudions des operateurs de type Kramers-Fokker-Planck dans la lim- 
ite semi-classique quand I'exposant du maxwellien associe est une fonction 
de Morse avec un nombre fini ng de minima locaux. Sous des hypotheses 
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supplementaires convenables, nous montrons que les premieres ng valeurs pro- 
pres sont reelles et exponentiellement petites et nous etablissons leur asymp- 
totique semi-classique complete. 
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1 Introduction 

In this article we shall continue the study that we started in [llj of the exponentially 
small eigenvalues of the Kramers- Fokker- Planck operator 

P = yhd,- V'{x) ■ H + |(?/ - hdy) ■ {y + hdy) on M^^ = Rf x MJ, 7 > 0, (1.1) 

and similar operators. Suppose that the potential V G C°^(]R'^; M) is a Morse func- 
tion on W^ such that 

(9°1/ = Oil), \a\ > 2, and \VV\ > l/C, for \x\ > C > 0. (1.2) 

Then (as we shall review) P is maximally accretive and has a unique closed 
extension from 5(M^") to an unbounded operator: L^ — )■ L^ that we shall also denote 
by P. The spectrum o"(P) is contained in the closed right half-plane. Assume for 
simplicity that 

V{x) -> +00, \x\ -> 00. (1.3) 

Then the Maxwellian e-(2^'/2+^(^))/'» belongs to the kernel of P, so G a{P). In p] 
the eigenvalues in any band < Kez < Ch were determined in the limit /i — )■ 
modulo 0{h°°), for every fixed C. They are of the form fih + o{h) where /z G C 
are the eigenvalues of the quadratic approximations of P with h = 1 at the various 
critical points of V. These /z values are explicitly known and belong to a cone 
|Im/i| < (9(Re/i). The o{h) terms have complete asymptotic expansions in powers 
of h (and we need fractional powers when certain multiplicities are present). 

Assume that V has no local minima, mi, 1112, ..., m„g and ni critical points of 
index 1 that we shall call saddle points. Then precisely uq of the above eigenvalues 
are o{h) (i.e. with /i = 0) and they are actually 0{h°°) (as can be understood intu- 
itively by using truncations of the Maxwellian near the local minima as exponentially 
accurate quasimodes). 

As we shall review below, it follows from the analysis in [11] that these no 
eigenvalues are actually exponentially small. In that paper we were able to establish 
the exponential decay rate and a full asymptotic expansion of the prefactor for the 
smallest non- vanishing eigenvalue when no = 2. 

In this paper we treat the case of general no and our results are similar to those 
for the Witten Laplacian by Bovier, Eckhoff, Gayrard, Klein [3], [1], obtained with 
probabilistic methods and the ones by Helffer, Klein, Nier [6], with simplifications 
by Nier [15j and le Peutrec [13] with full asymptotic expansions, based on the WKB- 
analysis and Agmon estimates in [S]. 

In order to state one of our results completely we describe first a geometric frame- 
work slightly generalizing the procedure in [6], [15]: Let mi be a global minimum of 
V and put E'mi = R*^- For cr ^ 1, V^^{] — 00, cr[) is a connected relatively compact 
open subset of M.'^. When we decrease a, V^^(] — 00, a[) remains connected until we 
reach a critical value o"2 where one of the following happens: 

a) o"2 = 0(mi) and V(] — 00, cr2[) is empty. The procedure then stops. 



b) V~^{] — oo, (T2Q is the finite union of several disjoint components, E1UE2U...UEN, 
N > 2, wliere tlie labelling is chosen so that mi G Ei. For k > 2, let nik G Ek be 
a minimum of V\^ and write £^^2 = -^2, •••, E^^ = Ejy. Notice that none of the 
closures of the Ej can be disjoint from the union of the other closures, and that the 
intersections of their boundaries are finite unions of saddle points. a2 is the common 
value of V at those saddle points. For 2 < k < N, we put a (Em,,) = 0^2- 

In case b) we pick (successively) each of Ej and consider Ej fl V^^{\ — 00, cr[) with 
a decreasing from o"2 until one of the scenarios a) or b) appears. In case a) we stop 
(with that component Ej) and in case b) (say for cr = era < (72) we get a finite union 
of connected components. Choose a global minimum for each of the new components 
except for the one which contains the already selected minimum rrij. We continue 
in this way until all the local minima have been recovered. Then for each local 
minimum we have an associated connected component E^ of V^~^(] — oo,cr(m)[). 
We put 

^(m) = a{m) — 0(m). 

See Section |4] for a more detailed description of the procedure. 
One of our main results is 

Theorem 1.1 Under the above assumptions, the hq exponentially small eigenvalues 
Hi, ...,/ino o'^e all real and can be labelled in such a way that 

uniformly when h -^ 0, and with the convention that S'(mi) = +00, fii = 0. 

This result is valid also for the Witten Laplacian —Ay = dydy, where dy : = 
0-v/hj^^^v/h ^^^ extends those of |1], [0], in the sense that no generic assumption is 
done on the separation of critical values. 

Under a generic assumption in the spirit of |1], |6], we have full asymptotic 
expansions for each of the eigenvalues /i^. 

Theorem 1.2 Suppose that V satisfies the same hypotheses as in the preceding 
theorem, and assume in addition the following : 

For every critical component E^ created in the above procedure, 

i) there is a unique point rrifc in Ek where inf^;^ V is attained, so Ek = E^^, 

a) except in the case when Ek = M.'^, there is a unique (saddle) point Sj = Sj(^k) in 
dEk which is also on the boundary of another component o/K~^(] — 00, a{Ek)[). 
In particular we always have N = 2 in case b) of the above procedure. 

Then we have the following asymptotic expansion 

Hk = hlk{h)e-^^'^^ with lk{h) = lk,o + lk,ih + ...,lk,o>0, (1.4) 

where 

, _ Ai(s,) /detr'(m,) 

Here — Ai(sj) is the unique negative eigenvalue of the block matrix ■ -.^„, ■. 
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See Theorem 5.10 and Section pi Under a weaker generic assumption (see The- 
orem 7.2) we get this for the lowest no n- vanishing eigenvalue /X2- This last result is 
also valid for the ordinary Witten Laplacian and seems to be new in that degree of 
generality. 

The proof is very much based on the analysis in the paper [11|, where we were 
able to treat the case of one or two minima. As there we shall also use in an essential 
way the supersymmetric structure of the Kramers-Fokker-Planck operator, due to 
Bismut pj and Tailleur, Tanase-Nicola, Kurchan [16], see also [T3] . 

What made it possible to go beyond the case of two wells was the observation 
that we have an additional (generalized PT) symmetry, namely if we conjugate 
the Kramers-Fokker-Planck operator P by the unitary and self-adjoint operator 
u{x,y) H- u{x,—y), then we get the adjoint P*. This extra symmetry makes it 
possible to introduce a Hermitian form on L^ for which P is formally self-adjoint 
and the restriction of this form to the spectral subspace corresponding to the Uq 
lowest eigenvalues is positive definite, hence an inner product. Consequently the 
restriction of P is self-adjoint and the no lowest eigenvalues are real. 

More generally, this extra symmetry entails that the supersymmetric approach 
followed in [llj reduces to a self-adjoint problem, very close to that of the Witten 
Laplacian, when we restrict the attention to the exponentially small eigenvalues. In 
particular, we can follow the work by Helffer-Klein-Nier [B] , who adapted probabilis- 
tic ideas in the case of the Witten Laplacian. Le Peutrec [H] simplified some of the 
linear algebra in [6], and we simplify that even further. 

The structure of the paper is as follows. In Section |2| we review some of the 
basic analysis of [11] , including a description of the supersymmetric formalism. The 
generalized PT symmetry is then explored in Section |3} first in the original Kramers- 
Fokker-Planck case, and then in the general supersymmetric setting. In particular, 
we show here that the first rio eigenvalues are real. In Section |4] we prepare for the 
multiple well analysis, by introducing a suitable labelling of the minima and the 
associated saddle point values, as well as a system of quasimodes adapted to the 
minima. Sections |5] and |6] are devoted to the analysis of the exponentially small 
eigenvalues in the generic case. The general case is studied in Section [7j The main 



results of this work are Theorem 5.10 Proposition 6/7, Theorem |7.1| and Theorem 

It is natural here to acknowledge a contribution of L. Boutet de Monvel. One of 
us had the priviledge to listen to his lectures at the Mittag Leffler institute in 1974 
about complexes of pseudodifferential operators and related operators with double 
characteristics, cf [2]. This type of problems continues to be important, and the 
present paper deals once more with that situation. 



2 Review of some results from 11 



The purpose of this section is to review some basic results from [TT] . We refer to that 
paper for proofs and more details and we give proofs only for slightly new variants. 
As in [11] some of the analysis works also on compact manifolds but currently we 



can go all the way only on M" and we restrict the attention to that case from the 
beginning. 

2.1 The general case 

On M" we consider a second order differential operator 

n 1 " 

P = ^ hD^^ o bj^k{x) o hD^^ + - ^{cj{x)hdcc^ + hd^^ o Cj{x)) + pq{x) 

j,k=i i=i (2.1) 

■' I OXj 

where the coefficients bj^k, Cj, po are assumed to be smooth and real, with bj^k = bkj- 
To P we associate the symbol in the semi-classical sense, 

Pix, = P2ix, + mix, + Poix), (2.2) 

n n 

P2{x,0 = Yl ^jAx)^j^k, Pi{x,0 = 'Ycj{x)^j, (2.3) 

j,k=l j=l 

SO that Pj{x, ^) is a real- valued polynomial in ^, homogeneous of degree j. (It is well- 
defined on T*]R" and coincides with the Weyl symbol mod 0{h?) locally uniformly.) 
We assume that 

P2(x,O>0, Po(a:)>0. (2.4) 

We impose the following growth conditions at infinity: 

d^bj^kix) = 0{l), \a\ > 0, (2.5) 

d^cjix) = 0{1), \a\ > 1, (2.6) 

d^poix) = 0(1), |a| > 2. (2.7) 

Put 

n 

iy{x,d^) = Ycj{x)d^.. (2.8) 

1 

Let f{t) e C°°([0,oo[; [0,3/2]) be an increasing function with f{t) = t on [0, 1], 
f{t) = 3/2 on [2,oo[, f{t) < t. Put /,(t) = e/(t/e), and introduce the time Tq 
average of /^ o pg along the integral curves of u, 



I rTo/2 
(/e o Po)to = 7f7 feOPo° exp{tu)dt. (2.9) 

-fO J-Tn/2 



By an averaging (realized by means of weak exponential weights) we showed in 
[TT] the following result: 



Proposition 2.1 Let P be of the form (2.1), where bj^k,Cj,po are smooth and real 



and satisfy (2.2)-(2.7). Define {fe opo)ro as in (2.9). Let Cq > be sufficiently large. 
Then for every C > 0, put e = Mh with M > sufficiently large. Then there exists 
C > such that 

Upo + hY^uW < CiUpo + h)-^P - z)u\\ + /^^lkll{(/,opo>.o<^})' (2-10) 

for u E S, Re z < Ch. 

Combining this witli a very simple and direct a priori estimate, we also deduced 
that 

\\B-^hDuf < \\{po + h)-^P - z)u\\\\{po + hY^uW +C\\h^uf. (2.11) 

Using elementary coercivity estimates and pseudodifferential machinery, we got 
indU, 

Proposition 2.2 

n{P- z) = L2, Rez<0. 

Here P : L"^ ^ L"^ denotes the graph closure of P : 5(]R") — )• 5(]R"). 

Corollary 2.3 The maximal closed extension Pmax of P (with domain given by 
{u G L^\ Pu G L^} coincides with the graph closure (the minimal closed extension) , 
already introduced. 

Thus P is maximally accretive. See [9], [7] for earlier and closely related results in 
this direction. 

We shall now discuss some weighted estimates for P—z, leading to simplifications 
and improvements in Section 3 of [IT] . These improvements will be used later on in 
this section. 

Let A = A(p) G C°°(M2"; ]o, +oo[) satisfy the bounds, 

X,X-' = Oi{pf% {p) = il + \p\Y^ (2.12) 

for some fixed A^o ^ and assume that 

A = Oi(A), (2.13) 

in the sense that 

a"A = 0(A(p)-l°l), a G N^"". (2.14) 



In the proof of Proposition 2.2[ given in [llj, we checked that 

A = A,,^ = {epf 
satisfies these assumptions uniformly for < e < 1, when A^ G M is fixed. The 



discussion there for such particular A = A^^at generalizes to A in (2.12), (2.13) and 
we have with A := Op^(A) (when h is sufficiently small): 



• The symbol of A'^ is equal to A'^ + Oi{h'^\-^{p)-'^). 

• The symbol of [-P, A] is of the form ^{p, A} + OoiJi^X), where we write a = 
Oo{m) if dp a = C(m) for every a G N^". 

• Here, 

{p, A} = {p2, A} + Oo(A) = -d^p2 ■ d^X + OoiX) = Oo(A(p)), 
where the term Oo{X{p)) is real- valued. 
Combining these facts with /i-pseudodifferential calculus, we see that 

• The symbol of [P, AJA^^ is equal to f %^ + Oo{h'^). 

Before continuing the main discussion, we shall give a simplification of the main 
step of the proof of Proposition |2.2 which is to establish: 



32; G C, Rez < 0, such that if -u G L^ 
and (P — z)u = 0, then u = 0. 



(2.15) 



To see this, we choose A = X^^n with < e ^ 1 and with A^ < —2 fixed. Then, 
= (A(P - z)u\Au) = (A(P - z)A-^Au\Au) = {{P - z - [P, A]A-^)Am|Am). 
Here we take the real part and use that Opi-Xj^^j^) is formally skew-adjoint, to get 
= ((P2 + Po - Rez + Op^iOoih^)))Au\Au) > (-Re z + 0(/i'))||Anf . 

Assuming that— Rez; ^ /i^, we conclude that An = and hence that u = 0. □ 

The improvement in comparison with [TT] is that we do not only consider Re ((P— 
z)u\u) but modify this expression with weight factors, while the weight factors in 
[TT] enter at a later stage. 

If we let u E L^ and drop the assumption that {P—z)u = 0, the same calculations 
and Cauchy-Schwarz give 

{-Re z + 0{h^))\\Au\\ < \\A{P-z)u\\, (2.16) 

for A = Ae^Tv with N < —2, and this remains true for more general A = Op/j(A) as 



in (2.12), (2.13), provided that A = Ci((p)" 



Proposition 2.4 Let A = Op(A) with X as in (2.12), (2.13). Then there exists a 
constant C > such that ifKez < —Ch"^, (P — z)u = v, u,v E L"^, Av G L"^, then 



Au G L^ and (2.16) holds 



Proof. Choose -N >No + 2. Then ( |2.16[ ) holds with A replaced by Ae,ArA: 

{-Rez + 0{h^))\\A,^NAu\\ < \\A,^nA{P - z)u\\, 



(2.17) 
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uniformly with respect to e. Letting e tend to 0, we get \\Au\\ < +00 together with 



(2.16). □ 



The proposition can be generahzed by letting u,v E S', with Aqu, Aqv G L"^, for 



some fixed Aq = Opi^{Xo) with Aq as in (2.12), (2.13) 



We shall next recall the dynamical assumptions introduced in [TT] . Consider the 
non-negative symbol 

p(x,0=Po{x) + ^^^. (2.18) 

Hypothesis 2.5 Assume 

The set {x G M"; po{x) = 0, u{x,d^) = 0} is finite = {x\...,x^}. (2.19) 

Let pj = {x^ , 0) and introduce the critical set 

C = {pi,...,p^}. (2.20) 

Notice that ^1,^0,^2,]? vanish to second order at each pj. For functions q on the 
cotangent space, we generalize the earlier definition of time Tq average and put 



1 rTo/2 
{q)To = 7^ q o exp{tHpJdt. 

J-O J-Tn/2 



''To/2 
-To/2 

We introduce the following dynamical conditions where To > is fixed: 
Hypothesis 2.6 

Near each pj we have {p)to ^ -^\p ~ Pj\^i (2.21) 

1 1 ~ 

In any set |x| < C, dist {p,C) > -, we have (p)to(p) > ^ , C{C) > 0. (2.22) 

^ C{C) 

V neighborhood U of tt^-C, 3^ > such that V x G M" \ t/, 

To To 11 (2.23) 

meas ({t G [-y , y]; po{exptu{x)) > -}) > -. 

We know that P has no spectrum in the open left half-plane. By using an 
elaborate method of micro local exponential (but bounded!) weights we showed in 
pT] the following result. 



Proposition 2.7 For every constant B > there is a constant D > such that P 
has no spectrum in 

{z eC;Rez< Bh,\lmz\> Dh} (2.24) 

when h > is small enough. Moreover \\{P — z)^^\\ = Osih'^) for z in the set 



Let pj G C and let Fp^ be the matrix of the hnearization of Hp at pj (the so 
called fundamental matrix of p at the doubly characteristic point pj). According 
to (2.21), the time average of the quadratic approximation of p at pj along the 
Hamilton flow of the quadratic approximation of pi at pj is elliptic and takes its 
values in a closed angle contained in the union of {0} and the open right half plane. 
We could therefore apply a classical result to see that the eigenvalues of Fp^ are of 
the form ±Aj^fc, 1 < fc < n, when repeated with their multiplicity, with ImXj^k > 0. 

Put 

^(a;, = -J5(a;, «0 =P2+Pi- Po- 

Let Fq, Fp be the fundamental matrices of q, p at one of the critical points pj G C. 

Since 

I, ,, . d . , . d 



Hq{x,0 = -{Pe{x,r]) 



dx 



P'x{x,v) 



drj' 



with 7] = iC,, 



Fq and ^Fp have the same eigenvalues; ±^Afc, k = l,...,n (j being fixed) where 
Re (|Afc) > 0. Now q is real-valued and we can apply the stable manifold theorem to 
see that the Hq-fiow has a stable outgoing manifold A+ passing through pj such that 
Tp.A^ is spanned by the generalized eigenvectors corresponding to +-Afc, k = 1 



,n. 



We also know that A+ is a Lagrangian manifold and that q vanishes on A+. 

Let A_ be the stable incoming ifg-invariant manifold such that Tp^A'^ is spanned 
by the generalized eigenvectors of Fg corresponding to — ^A^, 1 < /c < n. In [llj. 
Section 8, we established a transversality lemma for A± which together with a de- 
formation argument led to the following result. 

Proposition 2.8 We have 0'|(O) > 0, 0" (0) < 0. 

Here (f)± is the generating function for the Lagrangian manifold A±. 
Let 

t~r(p,p,) = ^^A,,fc. (2.25) 



In [TT] we also obtained the following two results. 



Theorem 2.9 We make the assumptions (2.1)-(2.7), (2.19), (2.21)-(2.23), and re- 



call the definition of C in (2.20). Let B > 0. Then there exists Hq > such that for 



< h < Hq, the spectrum of P in D{0,Bh) is discrete and the eigenvalues are of 
the form 

p,,k{h) ~ h{pj,kfl + h^^^'^'p^jAi + h'"''-''Pj,k,2 + ...), (2.26) 

where the Pj^k,o ore all the numbers in D{0, B) of the form 



f^j,kfi 



n 



e=i 



i^j,kAj,i + 2*^ (P' Pi)' ^^^^ ^iMA ^ N, 



(2.27) 



for some j G {1,...,A^}, A^ = #C. (Possibly after changing B, we may assume 
that \pj,kfi\ 7^ B, Wj,k.) Recall here that ±Xj/ are the eigenvalues of Fp at pj. 
This description also takes into account the multiplicities in the natural way. If the 
coefficients Vj^ki in (2.27) are unique, then Nj^k = 1 one? we have only integer powers 



of h in the asymptotic expansion (2.26). 
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Theorem 2.10 We make the same assum,ptions as in Theorem, 2^ For all B, 
C > there is a constant D > such that 

D h 

II {z-Py^\\<^, for zeD{0, Bh) with dist (z, a{P)) > -. (2.28) 

Still with j = jo fixed, let 

1 " 1~ 

f^o = -Yl ^^^^ + 9*^ (^' P^o)^ ^i^'^ (2-29) 



be a value as in (2.27) and assume that /io is simple in the sense that (i^i, ..., i^„) G N*^ 
is uniquely determined by /iq. In particular, every Xi for which z/£ 7^ is a simple 
eigenvalue of Fp. Applying a classical construction we got 

/i(/l) ~ /l(/iO + /iyUl + /iV2 + --) (2.30) 

with uniquely determined coefficients /ii, /i2, ••• and 

a{x; h) ~ ao(a;) + hai{x) + ... in C°°(neigh (x^«)), (2.31) 

where aj(a;) = 0{\x — x-'''|^'"~^-'^+), m = ^ z/^ and a^ has a non-vanishing Taylor 
polynomial of order m, such that 

(P - /i(/i))(a(x; /i)e-'^+(^)/^) = C(/i°°)e-'^+(^)/'^ (2.32) 

in a neighborhood of a;-"'. Actually any neighborhood fi CC M" will do, provided 
that 

1) 0+ is well-defined in a neighborhood of Vt. 

2) Hg^ ^OonU\{x^°}. 

3) fi is star-shaped with respect to the point x^° and the integral curves of the vector 

field Z/+ := (vr^)^ {Hg, ), where 7r^((x,0) = x. 

1A4- 



We also know that /^(/i) is equal mod 0{h°°) to the corresponding value in (2.26 ). 



If 7 C -D(0, B) is a closed /i- independent contour avoiding all the values /ij,fc,o in 



(|2.27|), and 

the corresponding spectral projection, then, 

||vr/.^(xae-'^+/") - X^e-'^+Z^IU^ = 0(/i-) (2.34) 

if X e C^{fi) is equal to one near x^°. It follows that xae~'^+'''^ is a linear combination 
of generalized eigenf unctions of P with eigenvalues inside h'-f up to an error 0{h°°) 
in L^-norm. 

We next review some exponential decay results from Section 9 in [11] and start 
with the case when C is reduced to a single point: 

C = {(0,0)}. (2.35) 

11 



Let = 4>+{x) G C°° (neigh (0; M)) be the function introduced following Propo- 
sition 2.8 so that A+ = A,^ is the stable outgoing manifold through (0, 0) for the 



if„-flow. Recall that by Proposition 12.8 



0"(O) > 0. 
Moreover, there exists G E C°° (neigh (0,]R"); M) such that 

{d^q){x, (j)'{x)) ■ dxG x x^, G{x) x x^. 



(2.36) 



(2.37) 



Let ^cir) = {x e neigh (0); G{x) < r} for < r < 1. Outside the set ilciCoe), 
we put 

^=0_e^(G), (2.38) 

where g = g{G) = InG for G > Gqb, so that g'{G) = 1/G. Then 

q{x, ?(x)) <-^, xe neigh (0, M) \ ^^(Coe), (2.39) 

if Co > is large enough. 

We extend the definition oi ip to a full neighborhood of x = 0, by putting 



giG) = In(Coe) 



1 



Goe 



(G-Coe), for 0<G<Coe. 



(2.40) 



Outside a small fixed neighborhood of we flatten out the weight. Let fs{t) 



(5/(|) be the function introduced before (2.9). For some small and fixed 6o > 0, we 
put 

i^ = f5oW = f5o{<P-^g{G)) (2.41) 

which is also well-defined as the constant 36q/2 for large x. Finally, by the same 



averaging procedure as in the proof of Proposition 2.1, we can add a term 0{e) 
to ip, supported away from a fixed neighborhood of 0, such that if ipe denotes the 
corresponding modification of ip, we have the apriori estimate, 

h\\e^^/%\\ + h^\BhD{e^^/''v)\\ < 0{l)\\e^'/\P - z)v\\ + Oih)\\e^^/%\\na{Coe), 

(2.42) 
uniformly, for |Rez| < Gh provided that e = Ah for A large enough depending on 
G. 



Now let fi(h) = ^i^k{h) be an eigenvalue of P as in (2.26), (2.30) and assume 



that /io is given by (2.29) and is simple, as explained after that equation. Then ii{h) 



is a simple eigenvalue of P and is the only eigenvalue in some disc D{fi{h), h/Go). 



Let uy\iKBix',h) be the approximate solution given in (2.31), (2.32) and let u = 
T^h-yix^EKw) be the corresponding exact eigenfunction, where 7 = dD{fi, ^). Using 



(2.42) we established the following result in |llj : 

Theorem 2.11 a) Outside any h-independent neighborhood of 0, we have 

M, B^^hDu = 0(e-^/(^'^)) 
in L'^-norm. 
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b) There exists a neighborhood Q of 0, where 



-4>+{x)/h 



u{x; h) = {a + r)e 

\\r\\L2^n), \\BhDr\\L2^n) = 0{h°°). 



(2.43) 



Remark 2.12. If we drop the assumption (2.35) and allow A^ — 1 more points 
p2, ■■■iPn in C, then Theorem 2.11 is still valid, provided that all Pj^kfl in 2.26 with 
j > 2 are different from the value /io, associated to pi = (0, 0). 



[X- 



We now drop the assumption (2.35) completely, so that C = {pi, ...,pAr}, pj 
■^ 0), pi = (0, 0). We then have the following extension of Theorem 



2.11 



Theorem 2.13 We make the assumptions of Theorem 2.11 , with the exception of 
(2.35). Let MwkbI^;) be as in that theorem and as there, we put u = TTh-yi^wKB) , with 



Co large enough in the definition of 7 above, so that iJ,{h) is the only asymptotic 
eigenvalue of P inside h^y that we can associate with the critical point pi = (0,0). 
Then u is not necessarily an eigenfunction of P (due to possible eigenvalues of P 
inside a disc D{p{h),o{h)), associated to other points in C), but the conclusions a) 
and b) of Theorem |2.11| remain valid. 



This result was not stated in [TT], but follows fairly directly from Theorem 2.11 
and the earher results by some easy and standard arguments: 



Let us first "eliminate" C\ {(0, 0)}, by introducing P :- 



P+E?<yx{ 



for 



some small and fixed a > 0. Here x is a standard cut-off function to a small 
neighborhood of 0. 



2.11 



Then P fulfills the assumption (2.35), so the conclusions a) and b) of Theorem 
apply tou = 7r/i^(xMwKB)5 where tt^^ is the spectral projection associated 



to P, h-f. 

The resolvent identity {z — P)^^ 
implies that 



[z - P)-i + {z- P)-\P - P){z - P) 



-1 



u = u + 



27ri 



[Z - P)-\P -P){Z- P)-\u^KBdz. 



(2.44) 



h'y 



Using the a priori estimate (2.42) 



2.10 



2.3 



P) ^X^WKB, together with Theorem 

(P — P)(2; — P)~^X^WKB is exponentially decaying in L^, uni 



with P replaced by P and v = {z 
and Corollary 



, we first see that 
brmly for z & hj. 

Using Theorem 2.10 and Corollary 2.3| once more, we next get that u : = 
{z — P)~^{P — P){z — P)~^x^WKB is uniformly exponentially decreasing in 
T>{P). In particular, ||m|| + ||P2/i_Dti|| = 0{e~^/^^) and we have the same 



estimate with u replaced by the integral in (2.44). 



Theorem 2.13 now follows by combining the above facts for the two terms in 
(|2l4l). 
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2.2 The supersymmetric case 

The Witten approach has been independently extended to the case of non-elhptic 
operators hke the Kramers-Fokker-Planck operator in [15] (in supersymmetric lan- 
guage) and in [1] (in terms of differential forms). See also [13]. We here follow the 
presentation of Section 10 and 11 of [llj and refer to that work for more details and 
proofs. 

We start by a quick review of that in the semiclassical case, then we establish 
some basic facts about the principal and subprincipal symbols, especially at the 
critical points of the given weight function. Finally we add some growth conditions 



and a dynamical condition, so that the results of Subsection 2.1 can be applied. 



2.2.1 Generalities 

Let 

A : (M")* -> W (2.45) 

be a linear and invertible map. Then we have the real nondegenerate bilinear form 

[u\v)a = {a''A{u)\v), u,ve A\Wy. (2.46) 

If a : a'=(M")* -^ A^(M")* is a linear map, we define the "adjoint" a^'* : A^(M")* -^ 
A*^(M")* by 

{au\v)A = {u\a^'*v)A- (2.47) 

(In the complexified case, we take the sesquilinear extension of {u\v)a and define 
a^'* the same way.) 

If a; is a one form, we have 



^■5 



(u'^y'* = {Aioy, (2.48) 

where ^ and ^ denote the usual operators of (left) exterior product and contraction. 
li u,v are smooth k forms with suppw fl supp w compact, we define 

{u\v)a = / {u{x)\v{x))Adx 



and denote by a^'* the formal adjoint of an operator a : C^{W; A'^T*W) — > 
©'(M"; A^T*M"). We can consider 

d^^ : Co~(M"; a'=T*M'^) -^ C^iR""; a''T*W), 

acting coefficient-wise, and a straightforward computation shows that 

(hd,^)^'* = -hd,^. (2.49) 

Let G C°^(]R";]R) and introduce the Witten (de Rham) complex 

d^ = e-'^/'^ ohdo e^l^ = hd+ {d(j)y : C^{W; a''T*W) -^ C^{R''; A^'+^T*M"), 

(2.50) 
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with d^ = 0. We have 

n 

d^ = J2{hd,, + a,,0) o dx^, (2.51) 

1 
where hd^^ + dx^(l) ^-cts coefficient-wise and commutes with dx^, so 

n 

<* = J2^-hd., + d,J) o AidxkY. (2.52) 

1 

The corresponding Witten-Hodge Laplacian is given by 

- A^ = d'^'*d^ + d^df*, (2.53) 

and we have 

-Aa = ^(-/i^xfe + d^^(j))Aj^kihd^^ + d^^(l)) 

+ J2 ^hd^A^(f) o dx^Aidxk^. 



(2.54) 



Now write 

A = B + C, B' = B, C* = -C. (2.55) 



Then (2.54) gives 

j,k 

+ ^{{9x^(t))Cj^khd^^ + hd:^^ o Cj-fc o {d^^cj))) (2.hQ) 

j,k 

+ J2 '^hd.^d.J o dx(;A{dxk)^ . 

j,k 

Note that the last term vanishes on 0-forms, i.e. on scalar functions. To recover the 
Kramers- Fokker-Planck operator (cf [T6]), replace n by 2n, 

1/0 V 



Then (|2.56|) gives for 0-forms: 

7 
2 



2 V-l 7 



^i°^ = I E(-H. + 5..0)(H. + dy,<P) + hH^, (2-57) 



where 

is the Hamilton field of with respect to the standard symplectic form ^ dyj A (ixj. 
If we choose 

(l)ix,y) = ^y' + Vix), (2.58) 



we get the Kramers- Fokker-Planck operator (1.1) 



- Af =yhd^- V'ix) ■ hdy + ^i-hdy + y) ■ (hdy + y). (2.59) 
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2.2.2 The principal symbol of the Hodge Laplacian 

The principal symbol of — A^ in the sense of /i-differential operators is scalar and 
given by 



j,k 

j,k j,k 

The corresponding real symbol q{x,^) = —p{x,i^) is given by 

j,k 

j,k j,k 



(2.60) 



f2.611 



It vanishes on the two Lagrangian manifolds A±^. 
We define 

^±=H,, . (2.62) 

Using xi, ...,Xn as coordinates on A±0, we get 

Z/+ = 2 X Aj^kdc^Jd,^ = 2A(0'(x)) ■ d, (2.63) 

j,k 

y- = -2j2Aj,kd,,<f)d,^ = -2A\<l)'{x)) ■ d,. (2.64) 

j,k 

Let Xq be a nondegenerate critical point of (j), so that A^ and A_^ intersect 
transversally at (a;o,0). The spectrum of the linearization Fg of Hg at (xo,0) is 
equal to the union of the spectra of the linearizations 

z/° = {2A(f)"{xo)x) ■ d^ and z/° = -{2A'(j)"{xo)x) ■ d^ (2.65) 

of u^ and z/_ respectively at xq- Thus we are interested in the eigenvalues of the 
matrices Acp", A"(j)" . Here, A^cj)" = 0"~ (A0")*0" has the same eigenvalues as A(f)", 
and similarly (j)"A, (f)"A^ are isospectral to A(f)". Thus 

The eigenvalues of Fg are given by ± 2A,;, 

^ ^ . (2.66 

where Ai, ..., A„ are the eigenvalues of A(j) . 



(Here the notation is different from the one used prior to Proposition 2.8 ) 
We assume from now on that 

B>0. (2.67) 

In [TT] we established the following result: 
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Lemma 2.14 Letn{x,dx) = Mx-dx he a real linear vector field onW^ . Letn± G N, 
n+ + n_ = n. Then the following two statements are equivalent: 

(A) M has n+ eigenvalues with real part > and ra_ eigenvalues with real part < 0. 

(B) There exists a quadratic form G : M" — )■ M of signature (n+,n_) and a constant 
C > 0, such that 

fi{x,dx){G) > hx\\ X e M". (2.68) 

Using this lemma we got the following proposition at the nondegenerate critical 
point Xq that we assume to be zero for notational reasons. Here p^ is the quadratic 
approximation of p at (0, 0). 

Proposition 2.15 a) Assume that the matrix Acj)" of v^ has m± eigenvalues with 
± real part > 0, m+ + m_ = n. Then there exists a real quadratic form Q{x,C,) on 
M^" such that 



Rep'{{x,0+^eHg{x,0) > ^\{x,0\', {^,0 e 



i)2n 



0<e<l. 



(2.69) 



b) Conversely, assume that there exists a quadratic form Q such that (2.69) holds. 
Then Acj)" has n± eigenvalues with ± real part > 0, where {n^,n_) is the signature 

of no). 



The condition (2.21) implies the existence of Q as in a) of the proposition. 



2.2.3 The subprincipal symbol 



We next look at the subprincipal term in (2.54), i.e. the last sum in that equation. 
As we saw in ^Jj, it can be rewritten as 



2hJ2{<P" o A'){dx,)% 



(2.70) 



Now we restrict the attention to a nondegenerate critical point xo of and 
we shall compute the subprincipal symbol of — A^ at the corresponding doubly 
characteristic point (xo,0). At that point (f)" o A* : T*JSJ^ — )■ T*^W^ is invariantly 



defined and it is easy to check that (|2.70|) is also invariantly defined: we get the same 

.UJ* 



quantity if we replace dx 



I;--; 



ClXny Ox\)--)(^Xni 



by ui,..,Un, UJ 



is any basis in the complexified cotangent space and w^, 
tangent vectors for the natural bilinear pairing. 



15 ■■i^ni 



, w; 



where wi, 



, LOrt 



is the dual basis of 



Assume that the equivalent conditions of Proposition 2.15 hold and denote the 



corresponding eigenvalues (that are also the eigenvalues of 0" o A*) by Ai, .., A„ with 
Re \j > for 1 < j < n+ and with Re Aj < for n_|_ + 1 < j < n = n+ + n_ . The 



eigenvalues of Fp are then ±2i\j (in view of (2.66) and the isospectrality of Fp and 



iFq reviewed prior to Proposition 2.8) , so 



trF„ 



71+ 


n 

-E2A. 


lmp>0 


n++l 



(2.71) 
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The subprincipal symbol of the first term in (2.54) (at {xq, 0)) is equal to 



1 " 

5^A,,fe-{-^a+5x.,0,^O+5.,0} = -Y,A,,,<Pl, = -tT{A<p") = -Y.Xr (2.72) 

j,k j,k 1 

Simplifying the last term in ( |2.54[ ), we get the full subprincipal symbol at (xq, 0): 

n 

1 j 

and hence 

1 ~ " 

-trF, + Sp = -2J2 A. +2Y,ii<l>"A'){dx,))% 

n++l j 

The eigenvalues of ^A'p" o A^){dxj)^dxj on the space of m- forms are easily 
calculated, if we replace dxi, .., dxn by a basis of eigenvectors wi, .., u;„ of 0" o A^, so 
that 

if o A'){ujj) = X.ujj, 

and dxj by the corresponding dual basis vectors u*. (Here we assume to start 
with that there are no Jordan blocks. This can be achieved by an arbitrarily small 
perturbation of A, and we can extend the end result of our calculation to the general 
case by continuity.) We get 

J2ir ° A'){dx,)% = J2 A.^>^ (2.73) 

j j 

A basis of eigenforms of this operator is given by Uj-^ A .. A Uj^, I < ji < J2 < ■■ < 
jm < n and the corresponding eigenvalues are Aj^ + .. + Xj^. 
Then the eigenvalues of 

-trFp + Sp, acting on m forms 
are 

n 

2(A,, + .. + A,„ - Y, \)^ 1 < Ji < •• < Jm < n. (2.74) 

n+ + l 

We conclude that if m ^ n_, then all the eigenvalues have a real part > 0. If 
m = n^, then precisely one eigenvalue is equal to 0, while the others have positive 
real part and the corresponding one dimensional kernel is spanned by c<;„++i A ... Aaj„. 

2.2.4 A symmetry for adjoints 

In this subsection we are concerned with symmetry relations for the A, * adjoints. 
If D : L'^{n; A^r*fi) -> L'^{n; A^T*n), then a simple calculation shows that 

D = (D^'Y''*- 
This can be applied to — A^ and we get 

(-A^)^*'* = -A^., (-A^t)^'* = -Aa. (2.75) 
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2.2.5 Quasimodes and spectral subspaces 



So far, we only developed the formal aspects of the supersymmetric approach. Now 
assume, 

Hypothesis 2.16 : M" — )■ M zs smooth and satisfies 

d:<f){x) = 0{l), d:{{Bd,<f),d^<f))) = 0{l), \a\>2. (2.76) 

Moreover, (p is a Morse function and there exists a constant C > such that 

\(j)'{x)\>l/C, \x\>C. (2.77) 



We also assume from now on that A in (2.55) satisfies 



B>0. 



(2.78) 



(0) 



Then P = — A)^ satisfies the assumptions (2.4)-(2.7) and Hypothesis 
fulfilled with 

C = {(x\0),...,(x^O)}, 



2.5 



IS 



where x^ denote the critical points of 0. Indeed, (2.60) can also be written 



and the vector field u becomes 

d 



(2.79) 



u = 2C0; 



dx 



This means that Proposition 2.2 and Corollary 2.3 apply to P = —A\^. Notice also 



FTT 



,(0) 



that the principal part of — A^ is scalar so that this operator differs from the tensor 
product of — A)^ with the identity in some C^^ up a zero order term which is 0{h) 
in C^. Hence the above mentioned results apply to —Ay as well, provided that 



we strengthen the assumption on z in Proposition 2.2 to Re 2; < —Ch, for C > 
suitably chosen. 



We now adopt the dynamical assumption. Hypothesis 2.6 Then Propositions 



2^ |2T| Ig^Sj Theorem ^ [210| [27L3] apply also. 
Proposition 2.17 For z G C, Rez < 0{h), v E S, we have 



'-A'r^' - z)-'d^v = d^{-A^^ - z)-\, 






{ih 



(^+1)> 



when z ^ cT(-Ay ) U (t(-A^ 



^A -z) d^ v = d^' (-AV 



z f, 



when z ^ a(— A 



^'-'^)Ua(-A^P] 



(2.80) 



(2.81) 
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Proof. By unique holomorphic continuation it suffices to establish these relations 
for —Rez :» h'^. 

For such values of z, we can apply Proposition 2.4 to see that A(— A^ — z)^^v E 



W_.A-K,^r2 T.f „. _ ^_aW._.A-1 



L^ with A = Op^((p)). In particular, ^^(-A^^-^)-^^ G L^. Let u = (-A 
so that u, d^u G L^. 

From (—Ay — z)m = f we get, using the intertwining relations; 



A 



z ^v. 



W ^\.. - i a(^+i) 



d^v = d^{-A\' - z)u = {-A\ 



z)aA,u. 



By the equality of the minimal and maximal closed extensions (Corollary 2.3), and 
the fact that d^v, d^u G L^, we get 

d^u = (-A|4+^^ - zy^d^v, 



and we get (2.80). The proof of (2.81) is similar. □ 

We shall now discuss the action of d^f, and d^'* on generalized eigenspaces of 

— A^. Let 7 C C be a simple closed contour such that no eigenvalues of the 

quadratic approximations of the operators — A^ at the critical points of for h = 1, 

i = 0, ...,n, belong to 7. 

Then h'y fl a{—A]^) = for < /i ^ 1, so we can introduce the spectral 

projections, 

nw = -L / (, + AiVrf., 

and their finite dimensional ranges, 

EW = nW(L2), 

Since S is dense in L^, we can replace L^ by S in the definition of E^^^ and Proposition 
2l7| tells us that H^^+^^d^ = d^U^^^ on S. Consequently, 

d^: eW^e(^+i). (2.82) 

In fact, if M G E^^\ we can write u = Il^^\u), u E S, and then d^u = d^Ii^^^u) = 

n(^+i)rf^M G E('+'\ 

Similarly, 

df* : ^(^+1) ^ E^'\ (2.83) 

In what follows, we shall mainly consider — A^^ and —A^. Let uq be the number 
of local minima, mi, ...,771^0 of 0. Then, if x ^ C^(M"') is a standard cut-off to a 
small neighborhood of 0, the functions 

/r = h-'^^'xi^ - m,)e-(<^-^('"^-))/^ k = 1, ...,no, (2.84) 

are quasimodes of —A)^ with eigenvalue in the sense that 

ii/riixi, (/ri/r) = for M A ._,. 

- AjjVr = Oie-'/^') in Ll ^ ' ^ 

Using also the calculation for the subprincipal symbol above, we get, (as in [TT] ) 
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Proposition 2.18 IfC > is large enough, then — Aj^ has precisely Uq eigenvalues 
in D{0,h/C) (for h small enough) when counted with their multiplicity. These 
eigenvalues are actually 0{h°°). 



Let E^^^ be the spectral subspace corresponding to the eigenvalues in the proposi- 
tion and let 11'^°^ denote the corresponding spectral projection. Then e^ := 11^^^ f I 
form a basis in E^^^ such that 

/r - ef\ B^'hDifi'^ - ef) = 0(e-i/^'^) in L\ (2.86) 



This follows from Theorem 2.13 and its proof. Later we shall choose the quasimodes 
/^ more carefully, using more refined cut-offs. 

Let Til be the number of saddle points, i.e. critical points of index 1. 

Proposition 2.19 IfC > zs large enough, then fork > sufficiently small, — AJ4 
has precisely ni eigenvalues in D{0,h/C) (counted with their multiplicity). The 
corresponding eigenvalues are 0{h°°). We denote by E^^^ C L^ the corresponding 
spectral subspace. 



The first part follows from Theorem 2.13| and the calculation of the subprincipal 



symbol (see \TT] for more details). That the corresponding eigenvalues are 0{h'^ 
and not just Oiji^) is also the consequence of a standard argument: 
Let 

ff^ = h-'^'^a, (x; /l)e-'^+(")/^ (2.87) 

where 

aj(x; h) ~ 0.^,0(3;) + haj^i + ... in C°°(neigh (sj,]R")), aj^^Sj) 7^ 0, 

and </)+ X Ix — s^p solves the eiconal equation q{x,(j)'^) = in a neighborhood of 
Sj. We assume that /• is a quasimode so that in the sense of formal asymptotic 
expansions, 

-A?/f^=/^.W/r' l^Ah) = Oih'). 
Using the intertwining relations 

aJ)<'* = <'*A«, A^^^d, = d,A^}\ (2.88) 

we see that d>'*fj and d^fj are quasimodes to — A^^ and — A^ respectively, with 



the same eigenvalue Hj = o{h). However, from Theorem 2.13 and the calculations 
of the subprincipal symbol above, we know that these two operators cannot have 
any non-trivial quasimodes at Sj with an eigenvalue o{h). Hence d.'*fj and d^fj 

vanish in the space of asymptotic WKB expressions, and using that — A^ = d^'*d^+ 

dtf,d^'* , we see that fj,j = 0{h°°), as claimed. 

Recall that A^^^ is the stable outgoing manifold through {sj, 0) for the Hg-6.ow 
and that (p'^Asj) > by Proposition 2.8 (Similarly we have a stable incoming 



manifold A^_^.) Recall also (see Proposition 2.15) that the linearization of Hg, 
at that point has n — 1 positive eigenvalues and 1 negative eigenvalue. Its matrix 
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is 2A(j)"{sj), according to (2.63). We let K^,K^ C A^ be the corresponding stable 
outgoing and incoming submanifolds of dimension k+ = n—1 and fc_ = 1 respectively 
and recall that K+ C A^,^, K_ C A^_ and — (f>{sj) — 4>±j vanishes to the second 
order on 7ir^{K±). It is also clear that A^, A^^^ intersect cleanly along K±, so we get 
(cf [m Section 11]): 

0+ - (0 - (j){sj)) X dist (x, 7r,(K+))', - 0(s,) - 0_ X dist (x, 7r,.(ir_))2. (2.89) 



We also know from Proposition 2.15 and the isospectrality recalled prior to (2.66) 
that 0"(sj)y4* (which is the linearization of —^Hg. , where we also notice that q 

vanishes on A±^) has precisely one negative eigenvalue, that the other eigenvalues 
have real parts > 0, and that ajfi{sj) 7^ is an eigenvector corresponding to the 
negative eigenvalue. Here we recall that (f)"{sj)A and (j)"{sj)A^ are isospectral, since 
{(l)"{sj)A)^'* = 0"(sj)A*. Also, if a*Q 7^ is an eigenvector corresponding to the 
negative eigenvalue of (f)"{sj)A, then by ([11, (11.27)]), 

(4ol«i.o(5,))A ^ 0. (2.90) 



Here we point out that the construction of the quasimode in (2.87) starts by 
choosing ajfi{sj) to be a non-vanishing eigenvector of (f)"{sj)A^, corresponding to the 
negative eigenvalue. 

3 Generalized PT symmetry and consequences 

3.1 First remarks on the KFP case 

In this subsection, as a preliminary and pedagogical step, we show the reality of 



the small eigenvalues of the KFP operator (1.1), n = 2d. This corresponds to the 
general supersymmetric case with 

<l>{x,y) = y^ + V{x), ^ = l[\ ])■ (3-1) 

Here, V^ is a smooth real-valued Morse function on M^ with V" G C^, \V'{x)\ > 1/C 
for |a;| > C. Let mi,...,m„Q be the (non-degenerate) local minima of V, so that 
mi = (mi, 0), .., m„Q = (rrino, 0) are the local minima of 0. In this case the spectral 
subspace £'*^°) is spanned by the functions e^ , k = 1, ..,no, where 

ef^ = /i-te-sW-'^)-'^(-^'0))xfc(x, y) + 0{h^) in L^. (3.2) 

where Xk G C^(]R"') was defined before. 

Let K : M" — )■ M" be given by k{x, y) = (x, —y). Put Uf^u = uo k, u E L^(]R"), so 
that Uk is unitary on L^(]R'^) and also self-adjoint and equal to its own inverse. We 
introduce the Hermitian form 

iu\v)^ = iU,u\v), u,veL^{R^). (3.3) 
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In the following our operators will be real so we can restrict the attention to real 
functions and differential forms. Notice that 

p* = U^^PU^. (3.4) 

Consequently, 

{Pu\v)^ = {U^Pu\v) = iP*U,u\v) = {U,u\Pv) = {u\Pv)^, 

so P is formally self-adjoint with respect to the Hermitian form {u\v)k.. 

Proposition 3.1 The restriction of {■\-)k. io E'^'^'' x E'^^^ is positive definite uniformly 
with respect to h, for h small enough. 

Proof. Since o k = 0, we see that there exists a^ > independent of /i, such that 

(e^lef ). = au + Oik), {e^^^lef}). = Oih^) for k ^ k'. 

Hence for u = ^"^ -UfcC^ , we get 

no 

(u\u 



. = J2^ak + Oih))\u,\^>\\uf/C, 



and the proof is complete. □ 

In conclusion, 

Proposition 3.2 The restriction of P : E^^^ — )■ E^'^'^ is self-adjoint with respect to 
the inner product {■\.)k- 

In particular, the no eigenvalues of P with real part < h/C are all real. Note that 
they also are 0{h°°). 

Remark 3.3. The Maxwellian e~^'^^^^~^'^ I'^^l^ is an even eigenfunction with respect 
to y, associated to the eigenvalue 0. Let us show that no other eigenfunction in E^^^ 
than multiples of the Mawxellian can be even. 

In fact, assume the contrary, f/^-u = -u 7^ 0, Pu = fiu, and apply f/^ to the 
eigenvalue equation. Using that U^P = P*Uk, we get P*u = fj^u. (Taking the 
differences of the two equations, we get (P — P*)u = so that the Hamilton field in 
P annihilates u). Taking the sum of the two equations, we get ^{P + P*)u = fiu, i.e. 
^{y — hdy) • {y + hdy)u = jiu (for almost all x). Since jj, is not among the eigenvalues 
of the harmonic oscillator part, we conclude that -u = and get a contradiction. 
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3.2 Extra symmetry and self-adjointness in the general case 

Now we consider the general supersymmetric case and adopt the assumptions of 



Subsection 2.2 Let k : M" — t- M" be hnear and satisfy 

K^ = 1. (3.5) 

After a hnear change of coordinates we may assume that k, = Hjjn-d © (— Hjjd), for 
some d e {1, 2, ...,n}. Also assume that 

(j) o K = (f). (3.6) 

Then k maps the critical set of (p into itself. We assume that 

K,{x) = X for all critical points of (p of index or 1. (3-7) 

On differential k-ioims we define U^ = k* as the pull-back in the usual sense. On 
0-forms, we get 

UkU = U O K. 

On differential 1-forms, we get Ui^u = k^{uj o k,)^ where u o k means the 1-form 
obtained from u by composing the coefficients with k. More generally, for /c-forms, 
we get 

Using the fact that pullback and exterior differentiation commute, together with the 



invar iance (3.6), we get 

U-'d^U^ = d^. (3.8) 

We next assume that 

kA = A^k\ (3.9) 

which can also be written as kA^ = Ak^, since k^ = 1. Equivalently, kA (or Ah}) is 
symmetric. In the case of KFP, we have 

■* = K-i \)'^='h^'-^'^'^={l -i) P-i«' 

and we see that 

kA= - { _, 

2 11-7 



Proposition 3.4 The bilinear form defined by 

{u\v)a,k ■= {U,,u\v)a 
is a Hermitian form on the space of square integrable k-forms. 
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Proof. We have 

iu\v)A,. = {iA''A)iA\K'))uoK\v) = 

{{a''{Ak'))u\vok) = {{a\kA'))u\v o k) = 

{iA^K){A''A')u\voK) = {{A''A')u\iA''K')voK) = 

{u\{a''A){a''k')vok) = {v\u)a,., 

where complex conjugate signs are absent since we restrict the attention to real 
forms. □ 

The same type of calculation shows that 

{u\U^v)a = {U,,u\v)a^- 
Proposition 3.5 We have 

Proof. Let m be a fc-form and v he a {k + l)-form and consider 

{u\U^d'^ '*v)a = (U^uld'^ '*v)a^ = {d^U^,u\v)A^- 

Similarly, 

{u\d^'*U^v)A = {d^u\U^v)A = {U^d^u\v)A^, 



and we get the proposition in view of ( |3.8 ). □ 

It follows that 

f/«A^t = A^?7,, (3.11) 

which implies that A^ is formally self-adjoint for the Hermitian form (|)a,k- I^ fact, 
{Aau\v)a,k = iU^,AAu\v)A = iAAtU^u\v)A = 

(f/KM|(AAt)^'*f)A = {Uf,u\AAv)A = {u\Aav)a,k, 



where we used (2.75). 
We also have 

{d^u\v)A,K = {Uf,d^u\v)A = id^U^u\v)A = iU,,u\df*v)A = {u\df*v)A,K, 

which shows that d,'* is the adjoint of d^ for our Hermitian product (.|.)a,k- 

3.3 Posit ivity questions 



First we shall prove that the eigenvalues are real. We know from (2.56) that 



Re ( -A^°^ ) > 0, 



so all eigenvalues have real part > 0. Proposition 3.1 and its proof carry over without 
any changes, so we know that — A^-^ : E^^'^ -^ i?*^°) is self-adjoint with respect to the 
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inner product (.|.)k = {\)a,k on E^^'. In particular, the eigenvalues of — Aj^ with 
real part in [0, h/C] are all real and = 0{h°°). 

We next consider — A . and a vector in E^^"* of the form, 



S- 



W = (/i-"/\.(x; h) + rj)e-'^+-^(^)/'^^j(x) + F^-e"'^"/'^, 



(3.12) 



where 5*0 is a positive constant, r^, B^^'^hDr.j are 0{h°°) in L^ and similarly for 



J' 



6'j G C^(neigh(sj)), 6j = 1 near Sj, and aj{x;h) ~ aj,o(a;) + haji{x) + 



where aj^ol-Sj) 7^ is an eigenvector corresponding to the negative eigenvalue of 
0"(sj)A*. The vect or e^ is obtained by the spectral projection of the truncation of 
a quasimode as in (2.87). As remarked at the end of Section |2l when constructing 
a quasimode for — A^, we consider a*Q 7^ 0, which is an eigenvector corresponding 
to the negative eigenvalue of (f)"{sj)A. Then 



djQ^Qj 



;oisj))A ¥" 0. 



(3.13) 



Notice that we can take ctjo ~ '^^^jfii^j)- 

Now fix j for a while, suppress this subscript from the notation and write Oq 
ajfi{sj). We shall study the sign of 

{k ao|ao)A, 
which we already know is real and non-vanishing. Recall the relations 

K*0"k = 0", 0" = (j)"{Sj), 

and that A = B + C with B = B^ >0,C = -C\ A bijective. Write 



(3.14) 

(3.15) 
(3.16) 



K 



1 

-1 



■.E, 



E^^E,® E^ 



(3.17) 



where Ej^ : 

iwy = E] 

k' - Id and El 



1" '^ is the kernel of k — Id and E^ ~ M'^ is the kernel of k + Id. Let 
E'l be the dual decomposition, so that E*, ~ (]R"~'^)* is the kernel of 



{Wy is the kernel of k* + Id. 



Then (3.15), (3.16) say that 



A 



$1 


^ 





$2 , 


3ii 


C12 


'-'12 


-B22 



E, 



E* 



E_ -^E* 



E*_ ^ E. 



E* 



E_ 



(3.18) 



(3.19) 



where Bu, B22 > 0. We now make a continuous deformation of A into the identity in 
such a way that the above properties of A are preserved, deform oq correspondingly 
so that it remains an eigenvector corresponding to the negative eigenvalue of cf)" A^. 



(Notice here that the dynamical assumption (2.21) is automatically satisfied once 
we make B positive definite.) Then the quantity (3.14) remains non-vanishing and 
hence of constant sign, so we have reduced the problem to that of studying the 
sign of (fi;*ao|ao) when Oq is a non- vanishing eigenvector of 0" corresponding to its 
negative eigenvalue. This gives 
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Proposition 3.6 There are only two possibilities: 



i) //$i is the component in (3.18) that has a negative eigenvalue, then the quantity 



(3.14) IS >0. 



a) //$2 is the component in (3.18) that has a negative eigenvalue, then the quantity 



(3.14) IS <0. 



A more invariant way of describing tlie two cases in the proposition, is to say that 
we are in the first case if 0" is positive as a quadratic form on E^ and that we are 
in the second case when 0" is positive on Ej^. In the case of KFP, we have 



V" 
1 



K 



1 

-1 ' ' 



(3.20) 



so we are in the first case and the quantity (3.14) is positive 



Recall the final observation in Subsection 3.2 Combining it with Proposition 



3.6 and the method of stationary phase as in the case of 0-forms, we get 



Prop osit ion 3.7 Assume that at every saddle point, we are in the case 1 of Propo- 
sition \3.6\ Then the restriction of {■\-)a,k to E^^^ x E^^^ is uniformly positive definite 
and if we equip E^^\ E^^^ with the scalar products {■\-)k o.nd {■\-)a,k, ihen the adjoint 
of d^ : E^^^ ^ E^^Us d^'* . 



Remark 3.8. Let us finally mention that our discussion applies to the case of the 
usual Witten complex. In that case, we let A = z : {W^)* — )■ M" be the applica- 
tion whose matrix is the identity in Euclidean coordinates. Naturally much of our 
machinery is redundant in that case, but it may be of interest to point out that 



our Theorems 'JA_, 1J1_ seem to be new in that degree of generality for the standard 
Witten Laplacian. 



4 Labelling and quasimodes for multiple wells 
4.1 Separating saddle points and critical components 



Let </) G C°°(]R";]R) be a Morse function satisfying (2.76) and (2.77). Also assume 
that 

^{x) — 7- +00, a; — 7- oo, (4.1) 



so that in view of (2.77), we have 0(x) > |a;|/C for |x| > C, where C > 0. In 



particular, e "^l G L . In [TT] we also treated the case when (4.1) does not hold 



and we hope to extend our present result to that case in the future. 

The function has finitely many critical points. The critical points of index 1 
will be called saddle points. In what follows, we shall be concerned with the saddle 
points and the local minima of 0. 
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Let s be a saddle point of and let B{s,r) = {a; G M"; |a; — s| < r}. Then for 
r > small enough, the set 

{x G B{s,r); (f){x) < 0(s)} 

has precisely 2 connected components, Cj{s,r), j = 1,2, with Cj{s,r) C Cj{s,r), if 
< r < r. 

Definition 4.1 W^e say that s G M" is a separating saddle point (ssp) if it is a saddle 
point and Ci(s,r) and 6*2(5, r) are contained in different connected components of 
the set {x G M"; 0(a;) < 0(s)}. Let SSP denote the set of ssp s. 

Notice that this definition depends on the global behavior of 0. It can be localized 
somewhat: 

Proposition 4.2 Let s eW^ he a saddle point and let o g]0(s), +oo[. Then s is 
a ssp if and only if Ci{s,r) and C2{s,r) are contained in different components of 
{x G F{s,a); (f){x) < 0(s)}, where F{s,a) denotes the connected component of s in 

r\]-oo,a[). 

Definition 4.3 A connected component E of the sublevel set 0"^(] — 00, a[) will be 
called a critical component (cc) if BE fl SSP 7^ or if E = M". 

We shall now describe a labelling system for the set LM of local minima of the 
function and a natural injective map from LM to the set CC of critical components, 
instrumental in constructing an appropriate system of quasimodes, adapted to the 
local minima. When doing so, it will turn out to be convenient to label the elements 
of both sets by 2-tuples of positive integers, N^. Let therefore mi 1 stand for a 
point of global minimum of 0, arbitrarily chosen, but kept fixed in the following 
discussion. Associated to mi^i we have the critical component Ei^i = M" and we let 
the associated saddle point value be ai = a{Ei^i) = +00. Let next 

(72 = sup0 (SSP) . 

Then the sublevel set C{a2) '■= {x G M", (pix) < (72} is the union of its finitely many 
connected components (all critical), of which there is precisely one containing the 
point rrii^i. The remaining connected components of the sublevel set C{a2) will be 
labelled as -E2,fc, 1 < ^ < -^2, -^2 > 1- Associated to each -E2,fc, we let m2,k be a 
point of global minimum of the restriction of to -E'2,fc, 1 < A; < A^2- 

Continuing the labelling procedure, we let a^ be the largest number of the form 
0(s), s G SSP, such that 0-3 < o"2- Decompose the sublevel set C^a^) into its 
connected components and perform the labelling as follows: We omit all those com- 
ponents that contain the already labelled minima rrii^i and m2,k, 1 < ^ < A^2- Some 
of these components may be non-critical. The remaining ones are critical and we 
label them as E^j, 1 < J < A'3, A^3 > 1. We then let rrisj be a point of global 
minimum of the restriction of to E^j, 1 < j < N^- 

We go on with this procedure, proceeding in the order dictated by the elements 
of the set 0(SSP), arranged in the decreasing order, until all local minima have been 
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a(l) 




mi 



Figure 1: Labelling in the generic case 



enumerated. In this way, we obtain the labelling of the set LM of local minima of 
by indices of the form k = {k'^, k'^'^) G N^, where k" is the index corresponding to the 
separating saddle point values a, arranged in the decreasing order, cti > 0-2 > o"3 > 
. . . > a^Q- It is then also clear that we get an injection from LM to CC, associating 
to each local minimum m^, the critical component E^^ containing rrifc. We equip the 
set of indices k with the lexicographical order. 



cr(l) = 00^ 




'7^1,1 m3^2 TO3,3 m2,i mi^i 7714,2 

Figure 2: Labelling in the general case 



Associated to the labelling procedure above is the connected graph, having a 
structure of a tree, or rather a root, obtained by letting a vary from +00 to —00, 
and representing the various components of C{a) as points that move vertically and 
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split at each separating saddle point level. 




(7(1) = 00 



a(2) 



a{3) 
a(4) 



Figure 3: Tree corresponding to the general case 



We also notice that our labelling procedure has the property that if -E^' and E^. 
are two critical components with Ek' C Ek, with a proper inclusion, then k' > k, 
and Ek can be reached from Ek' by following ascending links in the root. Given a 
critical component Ek, we let a{k) = (r{Ek) be the corresponding saddle point value. 
Here, as above, we adopt the convention that Ei^i = M" is a critical component, and 
that (t(1, 1) = +00. 



4.2 Cut-off functions and quasimodes 

In this section, we shall build quasimodes adapted to the local minima of and the 
labelling of the set LM, described in the previous subsection. When doing so, we 
shall first construct suitable cut-off functions, in terms of the corresponding critical 
components. 

Let E = E{a) 7^ R" be a critical component and let cr be the corresponding 
saddle point value. Let Xo G C^{B{0, 1), [0, 1]) be equal to 1 on B{0, ^). For 5 > 
small enough, put 

?(i) = ^{x) + V S%, (^) . (4.2) 



sessp, 

4>(s)><y 



6 
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For < t < -^, with C large enough independent of S, there is a unique connected 
component of 0^"'^(] — 00, c + t[) which is C6-close to E{a) in the sense that any of 
the two sets is contained in the algebraic sum of the other and -8(0, C6). Let Eq be 
this component when t = ^. Let / G C°^(]R; [0, 1]) be equal to 1 on ] — 00, a + ^] 
and have its support in ] — 00, a + -^[. The cut-off function associated to the critical 
component E = E{a) is defined as follows. 



XEix) = Ig fa;)/(0(a;)). 



(4.3) 
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Here 1^ is the characteristic function of the set £"0. 

Notice that outside the union of the balls B{s,6), for s G SSP, 0(s) > a, 
we have 0(x) > cr + ^, in the cut-off region supp Vx_b- We also notice that if 
B{s,6) n supp Vxb 7^ for some s G SSP, then dist {s,E) < 0{6). In particular, s 
can be a boundary point only of critical components E which contain E. 

We shall now describe the construction of quasimodes associated to the local 
minima nik, fc G N^, of 0. Here we assume that the minima have been labelled as 
described in the previous subsection. Let E^ be the corresponding critical component 
containing rrik, as described there. 

When mi 1 is a point of global minimum of 0, let Ei = M", and set 

When k 7^ (1, 1), we set, according to our labelling, 

for each of the local minima rrik, k G N^. Here XEk has been defined as in (4.3). 
Notice that then the quasimode fl is exponentially small near m^,^, as soon as 
k2 > ki. 

From the properties of the cut-off functions XEki 'we infer that our system of 
quasi-modes has the following two important properties: 

1) We have — Aj^ (/{^ ) = 0, and for the other indices fc, we obtain that 

Let (j{k) be the saddle point value associated to the critical component Ek and let 
Sk = cr{k) — (pimk). Then we get 

-^fUT)= E E lBM)[-Aj\x.J/^-?e-(^-^(-))/'^ 

-EfcCBeccsessPn^E (4.4) 

for some A^o > 0. 

2) We next claim that the quasimodes /^ enjoy the property of linear independence 
in L^, uniformly with respect to h. Indeed, notice first that ||/^, || x L Consider a 
linear combination 

where it is understood that the summation extends over all indices fc G N^, occurring 



in the labelling procedure of Subsection 4.1 For simplicity we now label them in 



lexicographical order, fki,fk2^ ■■■ fk„ , where uq = #LM, so that ki = (1,1). As 
already noticed, ii k ^ (1, 1), then /^ is exponentially small near mi^i = rukj, while 
fl has a substantial part of its L^ norm concentrated there: for every r > 0: 

\rk7rdx> ^ 



J Bi 



i/rprf^>^, c>Q. 
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Thus, if we multiply u by lB(mfc ,r)fk ^^^ integrate (assuming also that r > is 
sufficiently small), we get 



<-(0)m ^ Pfcil ^-At\^\„. |2 
\yu\-^B(mk^,r)J 

Hence 



IHiB(„..,,.)/r)l > ^-e"^(Ei«^'r'^- 



\ukA < C{\\u\\+e-cK\\u,\\), u. = M G M#^^ 

^^ \f'j^^\^dx = 0(e-i/(^'^)) for J > 2. Taking th 
we get 



Next consider the local minimum rrik^. We have j^,^ ^^ 1/^ \'^dx > 1/C, while 
^Bim,^,r) l/ifN^ = 0{e-'l^^^^) for J > 2. Taking the L^-product with li,(^,^,.)/£\ 



^ + 0(KJ) < C(hl| +e-c^||«.||), so Kl < C{\\u\\+e-M\u.\\). 

Continuing this procedure, we get Hm, || < CdlwH + e^^/*^'"^^||'U,||) and here the last 
term can be absorbed, so 



^Kr<o(i)ii^«,/, 



rr, 



which shows the uniform linear independence of the system of quasimodes /^ . 
Let e^ = n^''''(/^ ) where n*-°^ = 0{1) is the spectral projectio n of — A4 onto 



£;(o)_ Arguing as in Section 11 of PTI] (or as in the proof of Theorem 2.13) and using 



(4.4), we first see that 

ef = /W + O (/i-A^ie-(^^-°)/^) in L^ 

for some small fixed a > and A'^i > 0, and it is then clear that also the system 
(e^ ) is linearly independent in L^, uniformly with respect to h. 

In what follows, we shall need the quasimodes for the operator — A^ , associated 
to the saddle points of 0. The discussion here will be exactly the same as in the 
generic case, described in detail in the following section. 

5 Multiple well analysis in the generic case 

In this section we shall be concerned with the generic case and shall show that the 
analysis of the singular values of dcf, by Helffer, Klein and Nier in [6] in the Witten 
case can be applied here, thanks to the self-adjointness of — A*^°) in (.|.)a,k- We shall 
consider the case when e~'^/^ G L^(M") and assume from now on that we are in the 



case 1 of Proposition 3.6 for every saddle point s,. 



5.1 The critical points and quasimodes in the generic case 

Let be a Morse function satisfying the same assumptions as in the beginning of 
Section |4j The main hypothesis of this section is the following, 
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Hypothesis 5.1 For every critical component Ek = Em^ as in Subsection 4-i we 
assume that 

• 0|^ has a unique point of minimum (nik), 

• ifSSPnEk 7^ 0, there is a unique ssp s G SSPfl-Efc, such that sup (f){SSPnEk) = 
0(s) and in particular Ek fl (p^^i] — C)0, 4>{s)[) is the union of two distinct ccps. 



Combining tliis assumption with (3.7), we observe that K(mk) = ruk and i^{sj) = 
k{sj), for all 1 < A; < no, 1 < j < ui. 

In the following we let sq = oo. The general labelling procedure described in 
the previous section simplifies in the generic case, since here N^t = 1 for all k'^, and 
hence the elements of the sets LM and CC can be labelled by positive integers. We 
get the following result. 

Proposition 5.2 There exists an injective function 

{l,...,no} 3 k I — > j{k) e {0, ...,ni} with j {I) = 0, {sj(2),... ,Sj(no)} = SSP, 

and a family of connected sets Ek, for k G {1, ...,no}, such that the following pro- 
perties hold: 

i) We have Ei = M", and Ek is compact for k > 1. For every k > 2, Ek is the 
connected component containing rrik in 

{xGM"; 0(x) <0(sj(fc))}, 

and (f){mk) = min^;^. 0. 

a) If Sj G Ek and j = j{k') for some k,k' G {l,...,?7.o} and j G {l,...,ni}, then 
k' >k. 

Following [6], we shall now introduce suitable quasimodes, adapted to the local 



minima of and the simplified labelling, described in Proposition 5.2[ The construc- 
tion can be viewed as a special case of the general quasimode construction described 
in Section |4j In what follows, we let U stand for the set of the critical points of 0. 

Let ei > be such that the distance between critical points is larger than lOei, 
and such that for every critical point c and k G {1, ..., tt-q} we have either c E Ek ot 
dist(c,E^) > 10£i. 

When < Eq < El and Cq > 1 are to be defined later, we build a family of 
C^-cutoff functions Xk,ey k G {!,..., no}, < e < Eq, in the following way: first 
define xi,e = 1- For k > 2, we consider the open set Ek^e '■= Ek \ B{sj(k),£), and 
let Xk,e be a smooth function supported in Ek^e + B{0,e/Co) and equal to 1 in 
Ek^s + B{0, e/{2Cq)). Then for Eq small enough and Cq large enough, there exists C 
such that for all < £ < ^o; we have the following properties, 

(a) Xk,e is supported in Ek + 5(0,6:) and suppxfc,^ n {0 < 0(sj(fc))} C Ek- 
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(b) The distance from any local minimum and a separating saddle point c other than 
Sj(fc) to supp Vxfc,£ is bounded from below by 3ei. In addition, if c G suppxfc,£ 
then c & Ek. 

(c) There exists C^ > such that for all x G supp V(xfc,e)\-B(sj(fc),e), we have 

4>{sj{k)) + C^^ < (p{x) < (t){sj{k)) + Ce. 

(d) for any k' G {1, ..., uq}, if nik' G supp Xk,e then k' > k. In case k' j^ k we have 

0(mfc/) > 0(mfe) and (p{sj(_k')) < 4>{sj(_k))- 

(e) for any j' G {1, ...,ni}, if Sj G suppxfc,e then either j' does not belong to the 
image of the map j, or there exists k' > k such that 

nik' G suppxfc,e and f = j{k'). 

(f) Inside B{sj(^k),£) "we have: 

i) The distance from suppxfc,e H B{sj(^k),£) to the projection Tr^{K^) of the 
outgoing manifold K+ is bounded from below by a constant 6^ > 0, 

ii) for all x G B(sj(k),s) we have |0(x) — 0(sj(fc))| < Ce"^. 
The quasimodes associated to the minima rrik are introduced as follows, 

/W = /i-"/%(/i)e-^(*(^)-^("'=))xfc,.(a;), l<k<no, (5.1) 

where Cfc(/i) > is a normalization constant such that ||/^ \\f. = 1. Notice that 
Ck{h) ~ Ckfl + hck^i + . . ., with Ckfi 7^ 0. As before, to these quasimodes we associate 
their projections to E^^\ 

er = n(°)(/r). (5.2) 

In the same spirit we also define the quasimodes associated to the saddle points. 
For this we suppose that for all j G {1, ..., rii}, the cutoff function 9j is supported in 
B{sj, 2£:i) and equal to 1 in B{sj, ei). The corresponding quasimodes and projections 



are defined as in (2.87), (3.12). We call 



/W = (/i-"/\.(x, h) + r,-)e-'^+'^(")/'^%(x), (5.3) 

the quasimode, where aj(x, h) is the vector already introduced, so that ||/j \\a,k = 1- 
Again we denote by e\ the projection to E^^\ 



e 



(1) 



J 



n«(/f)). (5.4) 



It is of the form (3.12) for a new tighter cut-off 9j. 
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5.2 Estimates for the quasimodes 

We shall first work in the space E^^\ The following result says that, modulo expo- 
nentially small errors, the family (e^ ) forms an orthonormal basis in this space. 



Proposition 5.3 The space E^'^^ is spanned by the family (e^ )k=i,...,no defined in 



(5.1 5.2), and there exists a > 0, independent ofe, such that for all k, k' G {1, ..., no}, 



and < e < £0; we have 

{ef\ef), = 5,,y + 0{e~^'^). (5.5) 

Proof. We shall first repeat some arguments of Section 11 of pjj, essentially the 



same as in the proof of Theorem 2.13 Consider k G {1, ...,no}. Since 

/W = /^-"/4c,(/i)e-sW-)-*(-'=));^,_,(a;) (5.6) 

where Xk,e is the cutoff function defined above, we have 

P{fk^) = [P.XkA (/i-"/^c,(/i)e-sW-)-^(-^))^ 
= C>(/i-^Oe-(^*-^^)/'^), A^o > 0, 



(5.7) 



where Sk = <P{sj(k)) — 4>{i^k)- Here we used (d) and (f) of Subsection 5.1, where we 
gave estimates on the support of Vxk,e- Now we proved in Theorem 8.4 of [H], and 
recalled in Theorem 12.101 that 



n(o) 



^/(-pr'<i-o(i). ii(.-pni=o(i)...7. 



where 7 is an oriented circle of center and radius h/C with C large and fixed. This 
implies that 



e 



k 



n(0)^(0) ^ ^(0) ^ 0(/,-^ie-(5.-Ce)//.^ .^ ^2^ (5_ 



thanks to the following equalities: 






.(0)| .(0) 



where r^ is defined by (|5.7|). Now since fl is normalized in L^ in the sense that 



(/ri/r). = i>weget 

(erier)« = (/ri/f). + 0(/^-^^e-(^-^^)/'^) = 1 + 0(/i-^^e-(^-^^)/'^) 
and for all k, k' G {l,no}, we can write 



Let now k' 7^ k, and suppose that k' > k to fix the ideas. From Subsection 5.1 
we see that there are only three possible cases: 
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ruk' G suppxfc,£ (case (d) in Subsection 5.1). Then 

(/W|/W)^ = O (/,-^4g-Wm,,)-0(™.))//.) ^ 

since the support of Xk',e is included in the support of Xk,e- In that case we 
choose < a < (pijny) — (pirrik)- 

nik G suppxA:',e- This case is the same as above. 

The distance from nik' to suppxfc,e is larger than S^i (case (b) in Subsection 



5.1). Then by construction, the supports of Xk,e and Xk',£ are disjoint and 
hence (/ri/i°^). = 0. 

In all cases we can find a > independent of e and h such that 



Inserting this estimate in (5.10), and possibly shrinking Eq and a so that < a < 
(min(S'fc, Sk') - Ceq), we get. 



,(o)u(o) 



e';!,^), = 0{e-^'^), k'^k. 



The proof is complete. 



(5.11) 
D 



Proposition 5.4 The space E^^^ is spanned by the family (e- )j=i,...,ni defined in 

(5.3), (5.4), and there exists a' > (independent of e) such that for all j,j' G 
{l,...,ni}, 

{ef'>\e^^^)A,. = S,,r + 0{e'"'/'). (5.12) 

and aj{h) = aj^o + haj^i + .... 



Proof. The proof of this result is simpler than that of Proposition 5.3, thanks to 
the localization properties of the cut-off functions 6j. Recall that in (5.3) we defined 
the functions 

/(^) = (/i-"/4a,-(a;; h) + rj{x))e'^+'^^'=^/''ej{x), (5.13) 



(i)i 



,(1) u 



where aj{x,h) = Qj^q + haj^i + ... and ||/- \\a,k = 1- Again we denote by Cj its 
projection onto E^^\ Being in the case i) of Proposition 3.6, we can choose the 
coefficient a^ o such that 

iK^ajfl{sj)\ajfi{sj))A = 1. (5.14) 

It follows from [llj that 

p(i)(/«) = C(/i-^ie-"'/^), a' > 0, A^i > 0. 
Here we may also recall that 0+j(x) x |x — Sjp near Sj. Now using again Theorem 



8.4 of [TT] in the matrix case, or Theorem 2.131 we have 

n(i) 



^A.-p<')rd.=o(i). 
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where we recall that 7 is an oriented circle of center and radius h/C with C large 
and fixed. This implies that for some A^2 > 0, 

ef = ui')ff = ff + 0{h-''^e-^''^) in L', (5.15) 

Since the functions /• are L^-normalized (for the A, K-innei product), we get 
for all j, f e {l,...,ni}, 

{efVf^)A,. = (/ri/j'^)A,« + Oih-'^^e-'^'/^) = 5,,, + 0(/i-^^e-"'/^), 

where we also used that the supports of the functions 9j are disjoint, since the 
distance between different critical points is larger than lO^i. Replacing a' by a'/2 
completes the proof. □ 

5.3 Singular values 



Recall that in Section |2| as a consequence of Proposition 2.17 we showed that 



d(i, : E^'^^ —7- E^^\ Following [BJ, we shall now estimate the singular values of this 
map. The first step is the following proposition. 

Proposition 5.5 For all k G {1, ..., uq}, j G {1, ..., ^1} and < e < Eq, we have 
where bk{h) ~ bkfi + hbk,i + --- ■ Here h^fl is 7^ and will he determined in Proposition 



6.1. 



Proof. For the first result we notice that for all k G {1, ..., rio}, 

d^/r = {hS'\k.M) (/^-"/%(/^)e-^(^(^)-<^("^^») (5.17) 

is supported in supp Vxfe.e which is disjoint from the support of 6j for all j 7^ jik) 



according to point (c) in Subsection 5.1, As a consequence, the forms /j and d,^/^ 

also have disjoint support, and we get the first result, (/j \d^fl )a,k = 0. 

Now we consider the case when j = j{k). We shall use Gaussian type integrals 
in transversal directions, as in Section 11 of [H]. We first notice that the domain of 



integration in the inner product can be reduced to B(sj,e): from (5.17) and (5.3) 
we have 

(/j'V<^/r)A« = h'-'^/'cih) [ {{n\a,iKix),h))\S'^XkA^))^ + Oih^)) 

J B{sj,2ex) 

X dj o fi;(x)e-('^+'^°'^(^^+'^(^)-'^('"*^)/^t/x. 

(5.18) 
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Now we observe that for any x G B{sj, 2ei) \ B{sj, e) we have by (c) of Subsection 



5.1 that if X G supp Vxfc,e then 

(p{sj{k)) + C;:'^ < (p{x) < (l){sj(k)) + Ce. 



which imphes 

0+j O k{x) + 0(x) 



(mfc) > 



[X 



[rrik 



) = Sk + (p{x) - (t){sj(k)) >Sk + C^^. 



We can therefore restrict the domain of integration to the smaller ball B{sj,e), 
modulo an error term as follows, 



J B{sj,e) 



+ 0{h 



-N~{Sk+Cr^)/h\ 



(5.19) 



We shall need the following result. 
Lemma 5.6 In B{sj(^k),£), we have (p^jo k(^x) + {(j){x) —(f>{sj)) x dist (x, 7r^(-ft'_))^. 



Proof. Referring to the end of Section^ and especially (2.89), we first recall that 
0+(x) - (0(x) - (P{s,)) X dist (x, MK+)f, (5.20) 

and 



- 0_(x) + {(j){x) - (j){Sj)) X dist (X, TT^{K_)f 



(5.21) 



Now recall the definition of the adjoint operator — A^t = (— A^)"^ '*, from (2.75). 



Its principal symbol is p2 — ipi + Po = p{x, —0 ='■ Pix,0^ the corresponding real 
"g"-symbol is q{x,^) = g(x, — ^), and since this is the same type of operator, the 



geometric discussion above applies (notice in particular that Proposition 3.7 applies). 



We also notice for future reference that the hypotheses 2.5, |2.6| remain valid for 



— Ayit. Moreover, if we are in the case 1 of Proposition 3^ for — A^, then we are so 
also for — A^t . In particular we can introduce A^* , A^;^ as the stable outgoing and 
incoming iJ^-invariant Lagrangian manifolds through Sj and K^ C A^ as the stable 
outgoing and incoming manifolds for Hg. (noting that g = on A^^). Again we 

have dimKl = n — 1, dimK* = 1 and 



* ^■^2 



and 



+ 



[sj)) xdist(x,7r^(K:^)) 



(s,))xdist(x,7r,(K!))^ 



In view of the general relation 

J^iHg) = -Hg, where J : (x,0 ^ {x, -0, 



(5.22) 
(5.23) 
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we see that A,^;^ = J(A<^^), A,^. = J{A^_), or simply 



(5.24) 



From the relation (3.11) 



and Egorov's theorem, we see that p = po)C and hence g = go/C, where /C : {y, rf) i-)- 
{K{y), {K^)~^ri) is the natural canonical transformation which lifts n to the cotangent 
space. Since is invariant under composition with k, we have /C : A^ — )■ A^ and 
it follows that K^ = ]C{K±) and in particular ir^lK^) = k,{tt^{K±)). Combining 
([K22|), dpi, we get 



0+(x) + (0(x) - 0(sj)) X dist (x, 7r^(K:))2 X dist (x, fi;(7r^(i^_)))2, 

and replacing x by fi;(a;), we get the lemma. □ 

We summarize for further use some of the microlocal results obtained during the 
proof of the preceding lemma: 



Lemma 5.7 Let /C : (y, ■?]) 
(x, — ^). We have 



{n{y), (k*) ^7]) be the lifting of n and J : (x,^) 



A^*_ = J{A^+), A0. = J(A^_) and 0* 



y±, 



anc? a/so 



A^* = /C(A<^^), TT^i^K^) = K{n^{K±)), and (p*^ = (f)± o k 



(5.25) 



(5.26) 



i?n(i of the proof of Proposition 5.5 Now we are able to continue the computation 



in (5.19). Using Lemma 5.6, we get that the exponential term in the integral satisfies 



0+j O k(x) + (0(x) - 0(mfc)) - S'a: X dist (x, 7Tx{K.))^, 

on B{sj,e) fl suppVxfe.s- By evaluating the Gaussian integral in the directions 
transversal to n^^K^), we get 



(/j'^M,/r)A,« =hy%{h)e-'^/'^ + 0(/.-^e-(^^+^^-^)) 
=h'/%ih)e-'^/\ 



(5.27) 



It is also quite clear that bk is elliptic and a more detailed computation of bk.o will 
be given in Section |6j □ 

We are now able to compute the matrix of rf,^ from E^'^'^ to E^^^ with respect to 



,(0)^ 



,{1 



the bases (e^ ) and {e, ). This result is contained in the following proposition. 
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Proposition 5.8 There exists a" > such that if eo is small enough then for all 
k G {2, ..., rio}, j G {1, ..., ni} and < e < Eq, we have 



,(1) 



d^el 



(0)^ 



M,K 



O e 



-(^^+-")/'^) ./ j^m, 



and (ej.j^^lrf^ef )a„ 



h'/%{h)e-'''/\ 



(5.28) 



,(i)| 



M, 



where hk{h) ~ 6^0 + ^^fc,i + •■• is elliptic for k ^ 1, and (ej l^^e^ )yi^K = when 
k = l. ' ' 

Proof. Let k and j be integers as above. The case A; = 1 is clear since d^e[ = 0. 
Consider now the case when k > 2: using the self-adjointness and intertwining 
properties of the projectors 11*^'), / = 0, 1, which follow from Proposition 2.17 we 
check that 



.(i)i 






,. = (e«M.n(°)/rw = (e«|n«d,/r)^,. 



(n^^^efV^/i 



{0)n 
k M,K 



efM-> 



(o)> 



J \"'<PJk /^A 
'f-fl'^Mk')A,.. 



(5.29) 



(/«M,/r)A,.+(e«-#M./r 



According to Proposition |5.5[ the result of the proposition will follow if we show 
that for any k and j, we have 



(ei-'^-/f^M0/rk. = Ofe-(^^+"")/M. 



Let us therefore prove (5.30). By Cauchy-Schwarz 






(5.30) 



(5.31) 



Now a computation very similar to the one made in the proof of (5.7) and based on 
the support of Vxk e gives that 



d,fr = o {h-^o. 



No^-{Sk-Ce)/h 



) 



where C is the constant appearing in (c) and (f) of Subsection 5.1 On the other 



hand we have from (5.15) that 

ef - /f ) = 0(/i-^^e--'//^) in L^, 
where a' > is independent of e. Consequently, we can write that 

Taking ei small enough (e.g. such that Cei < a' /3) first and then posing a" = a'/2, 



gives the estimate (5.30) and the proof of Proposition 5.8 is complete. 



D 
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5.4 Main result in the generic case 

We can now give the complete asymptotics for the exponentially small eigenvalues 
in the generic case. Let us first prove the following result. 

Lemma 5.9 Consider the rii x Uq matrix 

R = i^j,k) ■= {{ej \d^e^ )A,K)jG{l,..,ni}, A:e{l,..,no}, 

where we recall that there exists a" > such that for all k G {2, ...,no}, and j G 
{l,...,ni}, 

rm,k = h^'Mh)e-'''^'\ r,,^ = 0(e-(^'=+"")/''), when j ^ j(^), 

while Tj^i = for all j . Set also ^j(i),i = 0. Then there exists rj > such that the 
singular values Vk{R) of R , enumerated in a a suitable order, satisfy 

Vk{R) = |r,(fc),,|(l + 0(6-^'/^^), 1 < A; < no. 

Proof. Since the first column of R consists of zeros, we know that vi =0 is a singular 
value of R and we can study the reduced matrix R' with entries r' ^ = Vj^^+i, where 
k > 1. We shall use that there is only one dominant term in each column of R'. 
Define the (no — 1) x (^o ~ 1) diagonal matrix D as follows, 

D = dia.g{h'/\+^{h)e-''+^/^ k = l,...,no - 1). 

Notice that D is invertible, thanks to the ellipticity of bk, and that VkiD) = \rj(^k),k\- 
Define the characteristic matrix of R' to have the entry 1 at each dominant term 
(columnwise), and elsewhere: 

U = (^i,i(fc+i))' 

where S is the Kronecker symbol. Then there is a constant rj > such that R'D~^ = 
U + Oie-"/''), 

R' = {U + 0{e-'^'^))D. (5.32) 

The Ky Fan inequalities therefore give 

yk{R')<{l + 0{e~^'^))uk{D). (5.33) 

To get the opposite estimate, notice that U is isometric, U*U = 1, and write 

U*R' = (1 + C(e-''/^))D, 

D = {1 + 0{e-"/''))U*R', 

to get 

UkiD)<il + Oie-'>/''))uk{R'). 

a 
One of the main results of this paper is the following theorem. 
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Theorem 5.10 In addition to the general assumptions, we adopt the Hypothesis 



5.1 and assume that we are in case i) of Proposition 3.6 at every separating saddle 



point. The exponentially small eigenvalues of P are real and given by 

lik = {h\h{h)\'' + 0(/^-))e-2S'=/^ A; = 1, ...,no, 

where Si = +oo (and ^i = 0), by convention. Here bk{h) ~ fo^^o + ^^fc,i + •••, where 
bkfl 7^ will be studied in more detail in Sectionl^ 



5.3 



and 



5.4 



the bases fe^ ) and (e^- ) are 



Proof. According to the propositions 

orthonormal up to exponentiaUy smaU errors in E^ and E''-^^ respectively, for the 
A, K scalar products. Let (e^ ) and (e^- ) be the corresponding orthonormalizations 
(obtained by taking square roots of the Gramians), which differ from the original 
bases by exponentially small recombinations. Then with respect to the ne w ba ses, 

and 



5.9 



the matrix of rf^ is i? = (1 + C»(e-"/'*))i?(l + 0{e-'^l^)) and from Lemma 
the Ky F an i nequlities (which will be used in a more essential way in the proof of 
Theorem 7.1) we see that the conclusion of that lemma is also valid for R. The 
matrix of the restriction of P to i?*^°^ with respect to the basis (e^ ) is i?*/2 and the 
theorem follows. □ 



6 Explicit computation of the leading tunneling 
coefficient 

The aim of this section to compute the dominant term in the amplitude of the ex- 
ponentially small eigenvalues. For this we shall estimate the coefficient b^ appearing 



in Theorem 5.10 We follow essentially the proof given in [8], but point out that 
here we are in a non-self adjoint situation. The main result is given in Proposition 
16.71 at the end of the section. 

6.1 Geometric preliminaries 

We recall some points concerning the study of the operator — A^. The principal 
symbol p(x,^) of — A^ does not depend on / and if q{x^^) = —p{x^i^) then from 
Subsection 12.2.21 we have 



q{x,0 = {A{^ + d4>{x)m - d<P{x))) = ((^ + d<P{x)m " d<P{x))) 



A, 



where we recall that A : (M")* — t- M*^, and where (■|-)a was introduced in (2.46). 
Near a critical point (s^, 0) of index 1, we also defined the two Lagrangian manifolds 
■^± = {{^lOi i = d(l)±{x)}, on which q is equal zero. This gives the following two 
eikonal equations in a neighborhood of Sj: 

= q{x,d(l)±j{x)) = {d(l)±j{x) + d(l){x)\d(j)±j{x) — d(j){x))A- (6.1) 
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We omit the index j when no confusion is possible, and define in a neighborhood of 
the saddle point Sj, 

^„ = -0_ + (0-0(s,)) ^^-"^ 

Notice that g- is the phase appearing in the Laplace integral defining hk in formula 
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5.19 up to a constant; indeed in a neighborhood of Sj we have according to Lemma 



o K + (0 - 0(mfc)) = (j)\ + {(j)- (pi-Sj)) + {(pi-Sj) - (pinik)) 

= -0_ + (0 - 0(s,-)) + (0(sj) - 0(mfc)) =g- + (0(s,) - 0(mfc)). (6.3) 



Using again Lemma 5.7 we also define 



t/; = 0; - (0 - 0(sj)) . . 

^:i = -0*_ + (0 - 0(s,)) = 0+ + (0 - <P{sj)) ^''- > 

and with these notations the eikonal equation for 0_ can be written near Sj in the 
following form : 

{dg\{x)\dg.{x))A = ^. (6.5) 

according to the following direct computation and (6.1[) 



{dgX{x)\dg-{x))A = ((i0+(x) - (i0(x)| - rf0_(a;) + (i0(x))A . . 

= ((i0_(x) + d(f){x)\d(f)^{x) — d(f){x))A = 0. 
Similarly using 0+, we get 

{dg*_{x)\dg4x))A = 0. (6.7) 

We now recall other properties of the functions g± and g'^, which are essentially 



reformulations of (2.89): In a neighborhood of Sj, we have 

g±{x) ^d{x,iT^{K±)f and g*^{x) ^ d{x,n^{Kl)Y (6.8) 

Recall some properties of the Hamilton fields i'± = Hq\^ near Sj that we 
identify with their x-space projections. We already used the fact that 

Vs e 7r,(ir±), z/+(x) = 2Ad(t){x) e T,(7r,(ir±), 

and that it vanishes at Sj. Similarly, 

u_[x) = -2A^d(f){x) 

so that z/_(x) = — z/^(x) and 

Vx G 7r,(ir;), z/_(x) = -2A'd(P{x) E T,{ti^{KI), (6.9) 

where the K^ are associated to the operator — Ay^t (see the end of Section |2]). 

When looking for an accurate expression of the eigenmodes of P = — A^ on / 
forms we need to know precisely the conjugate operator 

P^ = f,<P+/hp^-^+/h^ 

in a neighborhood of a saddle point Sj (recall that we omit the index j when no 
confusion is possible). 
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Lemma 6.1 In a neighborhood of Sj we have on l-forms 



=0 



P+ = hd{hd)'^'* + {hd)'^'*hd- {dg*_{x)\dg+{x))A (6.10) 

+ hCAdi4,+4,+) + h^Atd{4,-<t>+)^ 
where Cy denotes the Lie derivative in the v direction. 



tA,* I tA,* 



Proof. We first recall that P = d^pd,'* + d,'*d^ on forms (of arbitrary order) and 
we get 

P+ =e'f'+/\d^df* + <'*rf^)e-<^+/^ 

= {hd + d{(/) - 0+)^) {hd + d{<P + 0+)^)^ 



+ {hd + rf(0 + 0+)'') '* (/irf + rf(0 - 0+) 



-a[(p + (p+r'^'* 

yuu -r uytp — (p^j J . 

(6.11) 
For a one form u we have (cu^)"^'* = (Acu)^ so 

P+ = hd{hd)^'* + {hd)^^*hd 

+ rf(0 - 0+)^(M0 + 0+))J + (Arf(0 + 0+))Jrf(0 - 0+)^ 
+ hd{Ad{(j) + 0+))J + (Arf(0 + 0+))J/irf 

+ rf(0 - 0+)^(/irf)^'* + {hd)^'*d{(l) - 0+)^. (6.12) 

Using the formula (w, u) = z/Jw^ + w^z/J for a vector field z/ and a one- form w, we 
get 

rf(0 - 0+)^(Arf(0 + 0+))J + (Arf(0 + 0+))Jrf(0 - 0+)^ 

= (rf(0 + 0+)M(0 - 0+))a = -(rf(?l(x)|rf(?+(x))^ = 0. (6.13) 

Using also Cartan's formula CuCO = d^u^co) + u^dco with u = Ad^cj) + 0+) we can 
write 

hCAdin^P^) = hd{Ad{(t) + 0+))J + (Arf(0 + 0+))JM (6.14) 

and using also 

(zyJ)^'* = (A*-V)^, 

for a vector z/, we get 

= ((AV(0 - 0+))J)'''* (/.rf)^'* + {hd)^^* {{A'd{<P - 0+))J)^'* 

= d{(f) - 0+)^(/irf)^'* + {hdf'*d{(l) - 0+)^. (6.15) 
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Putting together (6.13 6.14 6.15) in (6.12), and using the fact that the eikonal equa- 
tion is satisfied by 0+, we get the lemma. □ 

It is convenient to introduce adapted coordinates in a neighborhood of the saddle 
point Sj. We begin with defining a C°° function yi by 

yi = ±\/^+M 

with the sign ± depending on a choice of a side of the (local) hypersurface 7ix{K+) = 
{gj^{x) = 0}. We impose for example the sign + on the side containing rrik for 
j = j{k) defined according to the injection defined in the preceding section. Of 
course we have 

yl = 9+ = <P+ - {<t>-<P{sj))- 



(6.16) 



We also introduce yl = yio k, so that according to Lemma (5.7) 

{yiy = {yi o ^? = 0; - (0 - <P{sj)) = g% 

and we notice that 

Ady{ E TTT.iK^) on 7r,.(K_). 

Now we consider the restriction of g^ to 7r^(i^+). Since 

5f_(.) X d{.,n^{K^)f near Sj, 

we can apply the Morse lemma to g^ on tt^IK^) and get smooth functions Z2, ■■■,Zn 
on Tix^K^) such that, 

9- = 4 + ■■■ + 4 on TTxiK+). 
Extend the functions z^ to a whole neighborhood of Sj using the Hamilton Jacobi 
equation 

^Ady-*Zk = 0, 

or equivalently 



{dyl,dzk)^ 



0. 



Since {dyl, dg^)^ = by (6.5), the function g_ does not depend on the yi coordinate 
and we have 

g^ = Z2 + ... + z'^ in a. whole neighboorhod of Sj. 

We introduce the variables zl = Zk o k, and we also have 

9- = {^2^ + ••• + i^nY in ^ neighboorhod of Sj. 



Note that according to (6.2) and (6.4) 



0+ = 1(9+ + 9-) = \{yl + i4? + - + {<?) , 

- 0(^.) = \{-9+ + 9-) = \ {-yl + {4f + ... + {zlf) 
Note that o k = gives 

-yl + {zlf + ... + {zlf = -{ylf + zl + ... + zl 

so that 

yl + zl + ... + zl = {ylf + {zlf + ... + {zlf. 



(6.17) 



(6.18) 



45 



6.2 A twisted Hodge operator 

We introduce now a special Hodge operator. We use the notation of [8]. We recall 
first the definition of the usual Hodge * operator, denoted * or X (and * or Xi when 
restricted on /-forms): For any real /-forms u and /i, 

{uj\^)dx = OJ A (*/i) = u A {Till), dx = dxi A ... A dxn, 

where (-I-) denotes the scalar product on /-forms inherited from the identification of 
TW^ and T*]R" via the standard euclidian metric. We denote by 

the corresponding application whose matrix is the identity in euclidian coordinates. 
This corresponds to a choice of an additional structure in our problem. 

In fact the right object is the following Hodge operator denoted *"^ or X^, and 
whose restriction on / forms will be denoted again by *^ or X/^: For any real /-forms 
u and /i, it is defined by 

{uj\iJ,)Adx = Lu A {*'^IJ,) = id A irfii). 

We give now some properties of the twisted Hodge operators *^ and *^ . 

Lemma 6.2 We have on l-forms 

i) *^ = *oa'(?^M*), u) (*^)-i = (-!)'("-') a' (^-M*)-^*, 

where * denotes X; in the first equality and Xn-i in the second one. In addition if u 
is an l-form and fi an {n — l)-form we have 

Hi) {*^uj) A (* V) = (det(z"^v4*))a; A fi. 

Proof. Let u and // be two /-forms. Then we can write 

u A *'^iJ, = {u\iJ,)Adx = (A^AuIjj,) dx = {u\ A^ A^jj,) dx 

= {uj\ a' (?"M*)/i) dx = UJA{*A^ («"^^*)/i) (6.19) 

This proves i). 

The invertibility of *^ and formula ii) are direct consequences of i) and the fact 
that X; oXn-i = (— 1)'("'~') for the usual Hodge operator. In order to prove iii) we 
recall first that 

{*u) A {*fi) = {*u\fi)dx = {fi\ * u)dx 

= /i A (* * w) = (-1)'("~'V Au = u Afi. 

Using property i) yields the result since 

{{A\t-^A)) A (A"-'(z-^A)) = A"(«-^A) = det(^-M)Id. 

The proof is complete. □ 

The second lemma is devoted to the commutation properties of *^ with the Lie 
derivatives. 
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Lemma 6.3 Let u be a smooth vector field. We have 
Proof. Consider the identity 



^A^A* = _£^ ^A ^^^ ^A*^A^ = _ ^A^ £^^ 



Cu{{u{x)\v{x)) Adx) = Cu{u A {*^v)). 

If M,f have compact support, then the left hand side is the Lie derivative of a 
compactly supported differential form and has therefore the integral zero. On the 
other hand it takes the form {Cyu\v)Adx— {u\Mv)dx for some first order differential 
operator M, so by integration, we recognize that M is the A* adjoint of Cy and the 
left hand side becomes 

{Cyu\v) Adx — {u\C^'*v)Adx. 

The right hand side is equal to 

The first terms in the two expressions coincide and hence the second terms also. 
This leads to 

-u A {*^C^'*v) = uAiCu *^ v) 

and the first identity in the lemma follows. 

In order to prove the second identity we shall first prove the general identity of 
independent interest: 

*^ *^' = *^%^ = (_i)'(«-0 det{i-^A') Id. (6.20) 

This follows from the computation 

det{i~^A^)uj A /i = (*^u;) A (*V) = (*^w|/i)A'^a; 

= (/i| *^ io)A^dx = 12 A (*^' *^ uj) = (-l)("-0'(*^' *^ uj) A /i. 

The second identity in the lemma now follows by applying *^ to the right and 



to the left in the first one and using (6.20). □ 

We also need a relation between *^ and n: 

Lemma 6.4 We have k**"^^ = det k*^ k*. 

Proof. The statement in the lemma is equivalent to the statement that 

K,*u A K* *^ fi = (det k)k*ui A *^k,*^ 

for all n — /-forms /i and all / forms u. Here the left hand side is equal to 

K*{u A * A*) = K*{{u\fi)A^dx) = (A'y4*a; o /tj/i o fi;)(det K)dx. 

The right hand side is equal to (det K,){K,*u\K,*fi)Adx so we only have to identify 
the scalar products in the two expressions: 

{k*uj\k*h)a = {a'' A a' K^tu o k\ a' K*/i o k) = {a\kAk^)uj o k\jj, o k). 

Here kAk^ = A^{k'^)^ = A^, so the two scalar products are equal. □ 
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6.3 Expressions for the quasimodes / 



(1) 



In this subsection we compute the leading amphtude of the quasimode /j in (5.3) 
on the manifolds 7Tx{K+) and n^^K^), associated to the saddle point Sj. 
As a warm up, we shall first show that 

*'^' {dz;A...Adz*^) = fdyi (6.21) 

in a neighborhood of Sj, where / is smooth and non- vanishing. In fact, composing 
the orthogonality relation {dyl\dzk)A = with k, it is not difficult to see that we 
also have 

{dy,\dzl)A^ = 0. (6.22) 

If a; is an n — 1 form, we have 

{uj\dz2 A ... A dz^)Atdx = u A *'^\dz2 A ... A dz*^). 



From (6.22) it follows that the left hand side vanishes as soon as u can be written 
as the exterior product of dyi and an n — 2 form, and then it is easy to see that 
(6.21) holds. f{sj) will be computed below. 

With these notations we have the following result: 



Proposition 6.5 In a neighborhood of Sj we have 

ff\x) = (/i-"/^a,(x, h) + 0{h^)) e 



-4>j^+{x)/h 



(6.23) 



where we recall that (f)j^+ = {yf + (^2) + ■•• + (^n) )/2 cii^d where 

aj{x, h) = ajfi{x) + haj^i{x) + ..., 

with Wx G 7r^(fC+), cijfi{x) = ajdyi 
and Wx G 7r^(-ft'* ), CLjflix) = aj * {dz2 A ... A (iz*). 
Here aj^ and aj^ are non-vanishing constants. 

Proof. We already know that ajQ{sj) has to belong to the kernel of |tr-Fp + Sp at 
{sj, 0) and it follows from the discussion at the end of Subsubsection 2.2.3 (or from 
[TT] ) that ajfi{sj) is an eigenvector in the negative eigenspace of (l)"{sj)A^. (See also 



(2.90).) However, dyi{sj) is such a vector since it is orthogonal to Ts^n^K^), the 
spectral subspace corresponding to the eigenvalues with real part > of A(f)"{sj) = 
{(/)" {s j) A^y . (Cf (2.66).) Thus we know from the start that ajfi{sj) is given by any 
of the two equivalent expressions in (6.23). 

Mimicking the proof given in [8], we use the expression of the conjugate operator 
P+ in Lemma 16.11 and have to solve 



P+ttj 



0, 



(6.24) 



in the sense of formal asymptotic expansions. We get the first transport equation 
Tajfi = 0, where 



rp dei ^ 



Ad{(j,+(l>+) 



+ c 



A,* 
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C 



Adg*_ 



c 



A,* 
A^dg+ 



where we used the definitions of g*_ and g+ in (6.4), (6.2). The transport equation 
for ajfl reads 



C 



Adg* 



-c 



A' 

AUg+ I "'J 



a,o = 0. 



(6.25) 



In order to prove the first identity in(6.23) we try to solve (6.25) along tTx{K^) 
with 

aj-o(x) = aj{z2, ..., Zn)dyi + 0{yi). (6.26) 

For general reasons, we already know that the vector field part of the transport 



operator is tangent to 7r^(-ft'+). Using (6.15), we get 

Here dg^ = on tTx{K^) so the first term vanishes there. Using the form of a^o 
above and the fact that dg+ = 2yidyi, we get {dg+)^ajfl = O^yf) and the second 
term also vanishes on iTx{K+). 

Thus we only have to solve along tTx{K+) the equation 



CAdg*_ajfl{x) = 0, 



(6.27) 



still with ttjfi of the form (6.26) and we have to check that aj is in fact constant on 
TTxiK^). From Cartan's formula for Lie derivatives, (6.27) reads 



= d [aj{z){Adgl_y dyi) + (^{Adg^ydaj) dyi — daj (Adg'ildyi) 

=0 

= d{aj {dg*_\dyi)A) + {dg*_\daj)Adyi 



(6.28) 



{dg*_\daj)Adyi = 2{d(f)\daj)Ady 



l5 



where we used that dg^ = 2d(j) on tTx{K^), and also the eikonal equation. Recalling 
that 2Ad(j) = z/+ G T{7i^{K^)) we get that on n^^K^) 



(6.27)^^VAd^5, = 



a, 



Cte, 



from the standard properties of non-degenerate vector fields. This proves the first 



assertion in (6.23). 



For the second one the tools are essentially the same. On t[x{K'^) we look for a 
solution of the type 

ajfl{x) = aj{yl) *^' {dz^ A ... A d^) + 0{z*). 

and we have to check that aj is in fact constant on tTxIK"^). 



Let us check that the first term in (6.25) vanishes on Trx{K^). dg*_ vanishes on 
'Kx{K*_), so 

^Adg*_ajfi = do {Adg*_yaj^o on 7r^(K* ), 
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and it suffices to check that 

d((Adg*_y *^' {dz;A...Adz*J 

i.e. that 



0, 



= di{Adg*_yfdy^) = d{fidg*_\dy^)A), 

and this is zero, since {dg^\dyi)A = as we have aheady observed and used. 
Thus the transport equation on 7r*(_ft'*) becomes 



■.A,* 



Ci+* iaM)d4^-^dz:) 



0. 



(6.29) 



Lemma 6.3 shows that the preceding equahty is equivalent to 

CA^dg+aj{yl){dz; A ... A dz*J = 0. (6.30) 

Exactly as in the preceding case and using again the eikonal equation, it reduces on 

71,{K*_) to 

2 {d(f), daj)^t dz2 A ... A dz'^ = 
and again this is satisfied if Sj = Cte: Indeed it also reads 

VAtd</,(aj(z/t)) =0, 

and 2A^d(j) = — z/_ G T{nx{K'i)) from (6.9). Again we used the standard properties 
of non-degenerate vector fields. The proof of Proposition |6.5 is complete. □ 

Now we evaluate the coefficients 5j and aj. For this we simply write that /• is 



[k, A)-normalized: 

= /z""/' f 9j o k{x)9j{x) (fi:*a,-(a;)|a,(a;))^e-'^+°"(")/^-'^+(")/'^da; 



where dj is the truncation function introduced in the preceding section 
Using (6.17) we check that 

1 



O K + 



2 {yl + {zlf + ... + {zlf) + I {{yir + z', + ... + 4) 

I {yi + zi + ... +zi) + l {{yir + iz;r + ... + iz:r) 



(6.31) 



*\2 



iyiy + i4y + ... + iz:) 



yl + Z2 + ... + zl, 



where we used (6.18) for the last two equalities. 

As for the amplitude in the integral, we use that at Sj we have two expressions 
for aj. To leading order w.r.t. h, we have: 



{K,*ajfi\ajfi)^dx = ajaj{dyl\ *^ {dz2 A ... A dz*^))Adx 

= ajaj{*^ {dz2 A ... A dz'^\dy\) A'^dx 

= ajCXj *^' {dz2 A ... A rf^:*) A *^'dyl 

= ±(det i~^A)ajaj {dyl Adz^A ... A dz^) 
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(6.32) 



where we used Lemma 6.2 and stopped trying to follow up the signs for simplicity. 



Putting (6.31) and (6.32) in the integral, using the change of variable n and applying 



the Laplace method, we get 

f(l)||2 
Ij \\k,A 



i=ii/.r^iiL 



so that 



±7r"/2(det i-^A)ajaj + 0{h). 



(6.33) 



We shall next compute f{sj) in (6.21). Equivalently, we shall compute g = g{sj 
in the relation 

*^ dyi = gdz2 A ... A dz*^. 

Indeed, if we apply *^* to this and use (6.20), we get 



det{i'^A^)dyi = g* dz^ A ... A dz* 



which is (6.21) with 



/ 



det{t-^A 



(6.34) 
(6.35) 
(6.36) 



(6.34) means that 

{u\dyi)Adxi A ... A dxn = ififw A dz2 A ... A dz^, 
for all 1 forms u (and with a sign independent of u). With u = dyl, we get 

{dyl\dyi)Adxi A ... A dXn = ^gdyl A dz2 A ... A dz^. (6.37) 



Now use (6.17) 



- <Pis,) = li-yl + iz;r + ... + {z:r) = ^i-iyiy + zl + ... + zl\ (6.38) 



where the last equality follows from 0ok = 0. On -K^iKS) we have Z2 
and hence 

Z/+ = 2Ad<\> = -2A{yldyl) = -2ylAdyl. 

Since i/+ is tangent to t^x^K^), we conclude that 

d 



z„ = 



so 



Adyl = b 



d 



dyi 



-2hyl— on 7r^(ii'_). 



(6.39) 
(6.40) 
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For notational reasons, we sometimes write zi, z\ instead of yi, y\. Recall that 



Zj {X} = Zj 



M^)) = y^ nj,kZk, 

fc=l 



where (nj^k) is the matrix of n with respect to the coordinates zi, ..., Zn- On nx{K^] 
we get yl = ni^iyi and (6.40) becomes 



d 
Z/+ = -2bKi^iyi-— on 7t^{K_ 
dyi 



(6.41) 



On the other hand, we know that 7r^(-ft'_) is an eigenspace of the linearization of 
Z/+ with associ ated eigenvalue — 2Ai < 0, writ ing A i = — Ai where Ai is the negative 
eigenvalue in (2.66), and a comparison with (6.41) shows that — 26/tii = — 2Ai, so 






(6.42) 



Combining this with (6.39), we get 



d A 

{dyl\dyi)A = {Adyl\dyi) = b {-—\dyi) = ^-. 

oyi f^iA 



(6.43) 



In Remark 6.6 we will give practically computable formulas for ki i. Also recall that 



the quantity (6.43) is > since we assume that we are in case i) of Proposition 3.6 



Hence ki^i > 0. From the last equality in (6.31) we also know that the matrix (Kj^k) 
is orthogonal and in particular that Ki^i < 1: 



< Ml < 1. 



(6.44) 



Inserting (6.43) in (6.37), we get 



Ai 



"«:i,idet^ 



(6.45) 



where we recall that we write zl for yl whenever convenient. From (6.31), we get 



1 

2' 



+ 6+0 k) 



fdz*Vdz* 
\ dx J dx 



so 



and (6.45) gives 



. 9z* , , , 1 
det^:— = ±(det -I 
ox 2 



O K 



\" U/2 



Al 



■«:i,i(deti(0+ + 0+o«:)^Ji/2- 



(6.46) 
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From (6.36) we now get at Sj: 



/ 



(det-( 



+ 0+°'^)L)'/'det(z-iA*). 



(6.47) 



Remark 6.6. Here is a direct way of getting ki^- Recall that ki^i is equal to yl/yi 
on 7ix{K^), so 

, 2 ivir 



'"''- yl 



on 7r^.(-ft'_). 



Using (6.38), we get [y]) = — 2(0 — 0(sj)) on nx{K^). Similarly, by (6.31), we have 



yt 



+ 0+ o K on Tix{K_), so 



K 



-(0-0(s,)) 



1,1 ~ 1, 



+ (b+O K) 



on 7r^(-ft'_). 



(6.48) 



It follows from the discussion that the right hand side of (6.48) is constant on 



TXx{,K^). Thus, in order to compute /tij up to the sign, it suffices to compute the 
eigendirection Tix{K-) and the Hessian of the positive definite solution 0+ of the 
Hamilton- Jacobi equation. 



Combining (6.21) and (6.23), we get 



(6.49) 



where / is given above, and using this and ( |6.47 ) in (6.33), we get 

1 = ±7r"/2 det{t-^A) det(z-^A*)^(det J(0+ + 0+ o fi:)L)'/'« •• (6.50) 



Here we notice that t : (M")* — )■ MJ^ is symmetric by its definition: {iuj\fi) = (w|yU.), 
a;,//e (M")*, so 

det^-M* = det(«~^A*)* = det A{i~y = detAt-^ = detz'M. 



Thus we get from (6.50) 



± 1 = 7r"/2(detz-iA)2(det J(0+ + 0+ o ^y^^y/^'^a'' 



Ai 



(6.51) 



where the sign to the left is the one that allows aj to be real. Assuming that ki i > 
(which can be arranged by the choice of sign oi yl), we get, 



a, = ±7r-4 (det2-^A)-i(det ^(0+ + 0+ o «:)'^J-4 



^ 

/^l,! 



(6.52) 



where we are free to choose the sign (implying a choice of sign for aj ) . 
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6.4 Computation of the coefficient 



Now we can compute the coefficient of the singular matrix defined in (5.27) according 



c(l) 



to (5.19) in the preceding section. First recall the expressions for the quasimodes 



c(0) 



/ and K in a neighborhood of the saddle point s.: : 



As for /■ we have the following definition from Proposition 6.5 in a neighbor- 



hood of Sj 



where we recall that 0j,+ = {yf + (^2)^ + ••• + {z^Y)/2. and that 
Va; G 7r^.(-ft'* ), 0^,0 (x) = cxj *^' [dz^ A ... A dz^). 



(6.53) 



Now we deal with the quasimode /^ where k = k{j) according to the injection 
defined in the preceding section. We have 

/(o) = /i-"/4c,(/i)e-sW-)-^(-'=))x,,,(x) 

where Xk,e is the cutoff function defined in the generic case and Ck{h) > a normal- 
ization constant. In the following it is convenient to pose 

Xk = v/det(0"(m,)) 

and a direct application of the Laplace method gives 



Ck = {l + 0{h))- 



n/4 



SO that 



f, 



(0) 



'Xh 



'l + Oih))e 



vr 



-{(j){x)-(l>{mk))/h 



(6.54) 



XkA^j 



(7r/i)«/4" 
We now choose the behavior of the cut-off Xk,e function near Sj. In addition to all 



the properties recalled in Subsection 5.1, we impose that in a small neighborhood V 
of 7r,(K_) n {y, > 0}, 

Xk,ei^) = XkAvi) and supp Xk,e n {yi < 0} = 0. 



As in Subsection 5.1 we impose that the function 9j is equal to 1 in a far larger ball 



of radius ^ ei ^ e. 



We can compute the tunneling coefficient defined in (5.19): 

= h^-'''^Ck{h) f Oj o k{x) {{K*aj) (x,/i)|dxfc,.(a;))^e-('^+-^°'^(")+'^(")-*('"'=))/'^cix 

+ 0(;i-^e-(^^+^^")/^) (6.55) 
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level line of 




We use the adapted coordinates. As already noticed the phase is up to a division 
by h equal to 

-(0+j o k{x) + 0(x) - 0(mfc)) = -(0+j o /t(x) + 0(x) - (t){Sj)) - Sk 

= -g-{x) - Sk 
= -{zl + ... + zl)-Sk. 
We therefore get 

Uj \d<j,fk )a,k 

= h}~^l^e-^'^'\k{h)( [e, o k{x) {{K*a,) (x, /i)|dxfc,.(x))^e-(^2'+-+^")/'^da; 



+ 0{h 



-N-(Sk+C7^)/h'. 



(6.56) 



where we recall that j = j{k). On 7r^(/C* ) we have by Proposition 6.5 
aj{x, h) = OLj *^* {dz\ A ... A dz*^) + Oih). 



Since z\ = z^ o k we can use Lemma 6.4 to get on 7Tx{K 



def 



Vx G 7ra,(_ft'_), Ci*jfl{x) = i^*Cijfl = (detfi;)Q?j * {dz2 A ... A dzn)- 
From the expression of Xk,e near Sj we also have in a neighborhood of tTx{K_ 

dxk,e{^) = x'k,e{yi)dyi, 



and putting these two expressions together with iii) of Lemma 6.2 , we get on 7rx{K_ 



{{K*aj){x,h),dxk,e{x))j^dx 

= ±ajXk,e{yi) {*^{dz2 A ... A dzn)\dyi)^dx + 0{h) 
= ±ajx'k,ehji) *^ {dz2 A ... A dz^) A *^dyi + 0{h) 

= ±aj{dett'^A)xk,e{yi)dyi A dz2 A ... A dzn + 0{h) (6.57) 
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Using the Laplace method in the z coordinates we therefore get from (6.56) and 



(6.57) 



= ±(det«~^A)/i 



\-nl2-Sklh 



ckmnhY^^''y% 



■K^K-fiV 



X'k,eiy^)dyi (6.58) 



= ±(det i-'A)7i^''-'y^h'/^e-^>'/^Ck{h)aj. 



modulo 0{h^^'^e •^fe/''). Combining this with (6.54), (6.52), we finally get 



Proposition 6.7 Let Sj be a saddle point corresponding to a minimum irik, with j = j{k). 
Then we have 



(/f^M./r; 



A,K 



± 



h 



TT 



(det0"(mfe))3 



fii,ij (det|(0+ + 0+o/t)"(sj-)) 



_e-5^/'>+0(/i3/2e- 



Sk/h\ 



where we recall that — Ai is the negative eigenvalue of A(f)"{sj) and that ^i^i is the re- 
striction of 

-2(0 -0(s,))' ' 



0+ + 0+ O K 

to T^x (yK^ ) (whose tangent space at Sj is the corresponding eigenspace). As a consequence 
the exponentially small eigenvalues of P = — A^ are real and asymptotically given by 



fik = h 



1 Ai 



(det0"(mfc))2 



TC 



«:i,i(deti| 



+ 



O K 



i"(^.)) 



r + Oih) e-2^'^■/^ fc = 1, 



,no 



where Si = oo (and fii = 0) by convention. 



To finish this section we shall make Proposition 6.7 even more explicit in the case 



of the Kramers- Fokker-Planck operator ( 1.1 ). Assume for simplicity that the saddle 
point Sj is placed at the origin x = 0, y = 0. Our calculations only concern the 
quadratic approximation, so we may assume most of the time that V is quadratic. 
After an orthogonal change of variables, we may assume that 



V{x) 



2 V~^ 



i^2 

V 



(6.59) 



where v,Vj are > 0. The variables can be completely separated and only the xi.yi 
variables are of interest, so in most of the following discussion we will assume that 
d = 1 with some remarks about the formulation of the corresponding results when 
d> 1. Thus we consider the one dimensional potential. 



V{x) 



-X 



(6.60) 
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^/-.2 , 2n 



^/-,2 „.2n 



where v > 0. Recall that 

P = ^iyi + vxri) + -^{y^ + ri^ 

q = -p(x, y; i^, ir]) = y^ + vxr] + ^{v^ - y^), (6.61; 

<P{x,y) = ^{y"^ -vx"^)- 

Let 0+(x, y) = ^{ax^ + 2bxy + cy'^) be the unique positive definite quadratic form 
which solves the eiconal equation 

q{x, y- d^(j)+, dy(j)+) = 0. (6.62) 

Expanding the left hand side as a quadratic form and equating the coefficients to 
zero, we get 

Ib^ + vb = 0, a + vc + 'jbc = 0, b+^c^ = ^. (6.63) 

Choosing 6 = as the solution of the first equation, leads to the two solutions 
±(f){x,y), neither of which is positive definite. Thus we have to choose the other 

solution, b = —2v/'j. Then the last equation in (6.63) leads to c = ±a/1 + ^ and 

the condition that 0+ is positive definite imposes the choice of the plus sign. Finally 

we determine a from the middle equation and get a = v^/l + ^. Thus the unique 

solution to our problem is 



V I Av o 2v 1 I Av 



which can be seen directly to be positive definite. (In higher dimesions, we get 

2\ 



(j)+\x, y) = (j)+\xi, yi) + ^ -(?/| + Vj 



+ 0+ O K = 


= "\ 


/. 4t; 2 


/. 4t; 2 



2 -J ■ 'j'^v 
■I 

differing from the expression for (f)^^\x, y) only in the variables Xi,yi, where we used 
the superscripts 1 and d to indicate the functions 0+ and in dimension 1 and d 
respectively. In this case, k : {x,y) i— )■ (x, —y), so 

/ /I -1. 

(6.65) 

Next we compute the negative eigenvalue — 2Ai of 2A(j)"{sj) and the correspond- 
ing eigenspace 7ix{K^). Since 

the eigenvalues are ^ ± a/(7/2)2 + t>, so 

Ai = ^( vW2F+^ - 7/2). (6.66) 
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the corresponding eigenspace 7ra;(i^_) is generated by the vector (1, a/(7/2)2 + v — ^) 
(In d dimensions, n^^K^) is generated by a vector {x,y) with {xi,yi) equal to the 
vector above and with the other components equal to 0.) 

nf ^ is the restriction of (— 20)/(0_|_ + 0_|_ o k) to 7Tx{K_) (where by our normal- 
ization, we have 4>{sj) = 0, Sj = (0,0)), so after some straight forward calculations, 
we get 

^ (6.67) 



Kii 



^/fT^' 



The formulae (6.66) and (6.67) remain valid in d dimensions, with —v denoting the 
negative eigenvalue of V"{sj) = V"{0). 

In the first formula in Proposition |6.7[ we write 



det 0"(mfc) 



det 



+ 



det0"(mfc) 
- det (p"{sj) 



X 



det (j)"{sj) 



det 



+ 0+ o k)"{sj 



(6.68) 



With V as in (6.59), the last factor reduces to the one-dimensional case and we get 

-det 



det 6" ( Si 



-V 
1 



1 



det|(0+ + 0+oft:)"(s,; 



det 



v^/l + f 





! + %■ 



(6.69) 



Combining (6.66)-(6.69) with Proposition 6.7, we get after some straight forward 
reductions. 






± 



h 



detV'inik) Vtl -^ , ,nfu"-^ - 



nj v-det1/"(sj 
so that the exponentially small eigenvalues are given by 

h f detV"{mk) 



h 



l^k 



n \-detV"{Sj) 



Sk 



Xie'~ + 0{h^)e~-. 



(6.70) 



7 Multiple well analysis in the general case 

In this section we return to the general case and show how the analysis of Sections 
|4] and |5] may be combined to obtain the other main result of this work. Theorem 
|7.1[ In Theorem 7^ we also give the full asymptotic expansion for the smallest non- 
vanishing eigenvalue under a generic assumptions which is weaker than the one in 



Theorem 5.10 The section is concluded by a brief discussion of an explicit example 
of a system with three minima and three saddle points, where the values of at 
the different minima and the saddle points may coincide, illustrating the transition 
from a degenerate case to a generic one. 



7.1 Statement of the main result and the matrix of df/, : 

In Section |4] we have constructed an injection from the set LM of local minima of 
to the set CC of critical components: rrik i— )■ E/., k G N^. Let a{k) = a{Ek) be 
the corresponding saddle point value of Ek-, with the convention that M" is a critical 
component and that cr(]R") = +00. It is of the form Ei^i where ttii^i is a point of 
global minimum of 0. For nik G LM, we put 

Sk = a{k) - 0(mfc) > 0, {Si = +00). (7.1) 

For simplicity, let us arrange the indices /c G N^ in a suitable order, to be chosen 
below, and write them as fci, k2, ■ ■ ■ kn^, with ki = (1, 1). 

Theorem 7.1 We label the local minima as above, so that Ei^i = M", and let ^2, •••, A^no 
denote the uq — I non-vanishing eigenvalues of —A]^ , which are o{h). For h small 
enough they are real and exponentially small. More precisely, with a suitable labelling of 
the eigenvalues, we have 

/if X /ie"25fe^/^, 2 < £ < no. (7.2) 

Next, we recall that in Section |4] we have constructed a basis for the subspace 
E(°), given by ef = n(°)(/f^), k e N^, and checked that the system (ef) is 
uniformly linearly independent in L^. Also, in Section [5| we have introduced a basis 
of one-forms (e^- ), 1 < j < ni, for the subspace Bp-\ associated to the saddle 
points, and this construction is applicable here as well. In particular. Proposition 
15.41 remains valid. 

Similarly to the generic case, we can analyze the matrix of d^ : E^^^ — )■ E^^^ with 
respect to the bases (e^ ) and (e^- ), and verify the following properties, 

Vi? = 0, /ci = (l,l), 

while for £ > 2, we have 

m 

i=i 
where for some fixed a > 0, 

1) r^.^ = hhj^e{h)e-^'=^/^, when Sj G SSFndEk,. Here bj^e{h) ~ E^=o^i,^^^ ^°/ ^ 0- 

2) r^-f = h^ (c)(e-('^(^^)-'^("'=^))/'^) + C(e~(^'=^+")/'^)), when Sj G SSP n dE and E ^ 
Ekf, is a critical component containing Ek^. 

3) rj^i = 0{h'i e~^^*'t^'^' '^) otherwise, i.e. when Sj is a saddle point that is not on 
the boundary of a critical component containing E^^. 



59 



The matrix of d^ : E^^'' — )■ E^^^ is of the form {OR') = R = (rj^i), where 
{j,i) G {1, ...,ni} X {1, ...,no} and indicates the first vanishing column. Let D be 
the (no — 1) X (no — 1) matrix diag (e~'^*<''')^=2,...,no; and write 

R' = h^RD, R = (Fj-^), (7.3) 

so that Tj^i is equal to bj^e, C(e~('^("^)"'^(™'=^)-'^'=^)/^) + C(e-"/'^), C(e-"/'^) respectively, 
in the three different cases above. Notice that in the second case, (p{sj) — (piruki) — 
Ski > t)ut if we want to allow cf) to vary nicely with parameters, we have no uniform 
bound from below by a positive constant. 

7.2 Singular values 

Clearly ||i?|| < 0{1) and we shall prove that R is injective with a uniformly bounded 
left inverse satisfying 

\\R-^\\<0{1). (7.4) 

Accepting this in this subsection, we shall now estimate the singular values of R. 
Using the standard notation, we let Si{R) > ... > Sn n-i{ R) be the singular values of 
i?, i.e. the eigenvalues of {R*R)^. Then in view of (7.4[), we have 



^ < s,{R) < C. (7.5) 



The Ky Fan inequalities and (7.3) tell us that 



s,{R')<hnR\MD), 

where we recall that ||i?|| = 0{1). Without loss of generality, we may assume that 
Sk2 ^ 'S'jtg > ... > Sk^ . Then s^{D) = e '="0+1-"' , and we get the upper bound, 

Su{R!) < O(l)/i5e"^'="o+i-/'^. (7.6) 



To get a lower bound we shall use (7.5) and write 

R = h-'^R!D-^ 
so that the Ky Fan inequalities give. 



^ < Sn,-i{R) < h-h,{R')sn,-AD-') 



Here SnQ-v{D ^) = e ^^0+^-'^' and we get 

Su{R!) > ^e-'''="o+l-/^ l<iy<no-l. 

In conclusion, 

s,iR') X /ile"^'="o+l-/^ (7.7) 
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We can now end the proof of Theorem 7.1 ,. We may assume that e[ , ..., Cni 
is an orthonormal basis in E^^^ for the scalar product {■\-)a,k- Let gl , ■..,gno be an 
orthonormal basis in E^'^^ for the scalar product {■\-)a,k on functions. Then, as we 
have seen, the eigenvalues fie are the squares of the singular values of the matrix 
R oi d(f) : E^'^^ — )■ E^^^ with respect to the two orthonormal bases. However the 
uniform linear independence of the basis e^ , ..., e^J means that R = RU where U 
and U~^ are uniformly bounded. By the Ky Fan inequalities we then know that the 
singular values of R and those of R are pairwise of the same order of magnitude. 



The theorem therefore follows from (7.7). □ 



7.3 Proof of (7.4) 



We define the indicator matrix R q of R by replacing Vj^e by in the third case in 



j,e) l<J<ni , 
2<l<nQ 



the description of R above (after (7.3)). Then 

Ro = {r 

where 

1) ^j,i = ^j/ = bj/ when Sj G SSP fl dE^^, 

2) rj,e = Tj^e = C'(e-('^(«^)-<^('"fc.)-^fe.)/^ + e'"/'^) when Sj E SSP n dE and Ek, C E e 
CC, Ek^ ^ E, 

3) Vj^e = otherwise. 

Then R = Rq + 0{e~°'^^) and in order to prove (7.4) it suffices to prove an a 
priori estimate 

||m|| <C»(l)||i?oM||, MeC"«-\ (7.8) 

since we then get an analogous estimate for R. 



We shall prove (7.8) by estimating the components of u one after the other in a 
suitable order and thus eliminating them successively. 

If SSP is empty or equivalently, if irii^i is the only local minimum of 0, then there 
is nothing to prove. 

If SSP 7^ 0, let o"2 be the smallest value in 0(SSP). Here we are not using the 
labelling introduced in Section [4| The components of 0^^(] — oo,cr2[) can be split 
into the critical components, labelled as Ek-, with k = {k'^^, k'^'^), 1 < k'^'^ < Nk'^2, and 
the others, say, Fi, ...,Fm- For notational convenience, we shall relabel the critical 

components E/^ a.s Ei, . . . Ej^. Here each Ej intersects at least one other, E^, while 
each Fi is disjoint from the closure of the other components. Thus the union of 
the Ej can be grouped into say, Ei U ... U Ee-^, -E^i+i U ... U Ee-^+e^,... such that 
the corresponding unions of the closures are connected and mutually disjoint. Here 
£i,£2,-- are all > 2. 

In what follows, it will be useful to change the labelling of the local minima 
whenever we find it convenient and here even the global minimum rrii^i may have 
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to change its notation. Consider one of the groups introduced above, say Ei, ..., Ee^ 
and write i instead of ii for simphcity. Let ruk be the unique local minimum of 0|^ 
(recalling that (T2 = inf 0(SSP)). Recall that in the construction of the quasimodes, 
we have assigned a critical component Ek to each local minimum and since (T2 is 
minimal, we have Ek D Ek- The same construction also shows that Ek 7^ Ek for 
precisely one value of k, say for k = 1. Notice that the corresponding local minimum 
mi is also a global minimum of the restriction of to £"1 U ... U Eg. 

Since EiU ...U Ei is connected, we know that EidEk ^ ^ for at least one value 
oi k y^ 1, say k = 2. This intersection is a finite set of ssps. Choose one of them and 
denote it by S2. We shall now estimate U2. Notice that r2,2 = ^2,2 is elliptic, while 
^2,1 = 0. Since S2 cannot be in the boundary of Ek for any A; > 3, we know (using 
also the minimality of (T2) that r^2,k = for all other k associated to any other local 
minimum different from mi and m,2. It follows that 

(^o«)(2)= 62,2^(2), 

and by the ellipticity of 62,2, we deduce the a priori estimate 

\ui2)\<0{l)\\Rou\\. 

If £ > 3, we may assume, after relabelling of m^, ...,m,e, that E2r\ E^ ^ 0. Let 
S3 be a ssp in this intersection. Then (using again the minimality of (12) we have 

(^om)(3)= 63,2^(2) + 63,3^(3), 

so, using the ellipticity of 63^3, we get 

|m(3)| < Oil)i\\Rou\\ + |n(2)|) < 0{l)\\Rou\\. 

Continuing this way, we get 

\u{2)\ + ... + \u{i)\ <C(l)||^o^||. 



Thus, in order to prove (7.4), it suffices to do so when u{2) = ... = u{i) = 0. 



Indeed, for a general u G C"° ^, write 

e 
I +'^u{k)ek, u{k) =0, foT2<k< 



u = u 

k=2 



where e^ is the canonical basis vector of index /c, and if we assume that (7.4) holds 
when u{2) = ... = u{€) = 0, then 

lln'll < Oil)\\Rou'\\ < Oil)i\\Rou\\ + 0(l)(|n(2)| + ... + |m(£)|)) < O(l)||^o^||. 

In other words, we have eliminated the minima 7712, ...,m,£ from the discussion 
and only m,i survives among the m,i,m2, ...,m,i. We do the same for the other 
groups -E^i+i U ... U Ei-^+i^, ..., (if there are more than one) so that for each group, 
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we eliminate all the local minima except one which is also a global minimum in the 
corresponding union. 

If (72 is the only element of 0(SSP), then we only have one group as above and 
no Fj. We have then eliminated all the u{k) and the proof is complete. 

If 0(SSP) contains more elements, let 0-3 > a2 be the smallest one. Again (after 
a change of notation and forgetting about the earlier Ek and Fy) we have 

rW - 00, (T3D = ^1 U ... U E;v U ^1 U - U Fm, 

where Ey are the ccs and F^ are the other components. Each of the components may 
have several local minima but we have already eliminated all but one, say m which is 
also a global minimum for that component, and Ey C Em and F^ C Em respectively. 
This implies that we can repeat the procedure of elimination precisely as we did at 
the level a2- After finitely many steps all the local minima are eliminated and we 



get (7.4) 



7.4 Full asymptotics for the smallest non- vanishing eigen- 
value 



We return to the general situation in Subsection 7.1 and label the minima tti^^, 
mfcj • • • if^kn so that irik^ = rrti^i and Sk2 is maximal among all the Sk-, j = 
2, 3, ..., no. We assume that there is a gap between Sk2 and the other Sk : 

Sk2 > maxS'fc, =: S'. (7.9) 

i>3 ' 

Also assume that 

dEk2 contains precisely one ssp. (7-10) 



Theorem 7.2 We assume (7.9), (7.10), in addition to the general assumption of Theorem 



7.1 Then the smallest non-vanishing eigenvalue ^2 is given by 

fi2 = ih\b2{h)\' + 0{h^))e-''^/', 



where we write S2 instead of Sk^far short. Here 62 is as in Theorem 5.10 

As before, we may assume that the basis {e, ) in E^^^ is orthonormal, while 

we have uniform linear independence for the basis (e^ ) in E^^\ We may assume 
however that 

l|e£^IU = l. (7.11) 



Define R,R',R,D as in and around (7.3). Then 



b2{h) Oie-"^!^) \ ( e-^2//j g 



««=(o(^/') V j( b)' <"2' 

where D = diag ie^^^^^)3<j<no and we write Sj instead of Sk for simplicity. 
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Here, we know that 62 is elliptic of order and that 
bounded left inverse. Consequently, 

1 



62 

W 



0{l) 



<W*W <0{l). 



has a uniformly 



(7.13) 



Theorem 17.21 follows from: 



Proposition 7.3 The smallest singular value ofRD : 



B 710-1 



-)■ 



P"l 



is equal to \h2\e 



-Si/h 



+ 



0{e (•S'2+a)/^) j'of some a > 0, independent ofh, when IR"^ is equipped with the standard 
Hilbert norm and M"°~^ is equipped with the norm [■], defined by 



\x\ 



E( (0)| (feJN 
XjXkyCj |e^ )k- 

2<j,k<no 



We know that [■] is uniformly equivalent to the standard norm and that [62] = 1, 
where 62 = (1, 0, ..., 0) e ]R"°~^ The square of the singular value is equal to 



inf \\RDxf. 

[x]=l 



(7.14) 



Eliminating the terms 0{e "'^^) in (7.12) will modif y the quantity (7.14) by a 



factor 1 + 0{e °'f^) and if we do that elimination, then (7.14) becomes 



inf \h2\^e-^''^"'xl + \\WDx'\\\ 

[x\=l 



X 



[X2, X 



(7.15) 



Taking X2 = 1, x' = 0, we see that this infimum is < |62pe~^'^^'''^, and we shall 
see that it is very close to this value. Let x be a point ( with \ x\ = 1) where the 
infimum is attained. Then ||iyi)x'|| = 0{e~^^^'^) and from (7.13) and the fact that 
3-^ = 0{e^'/''), we conclude that ||x'|| = 0(e-}^^-^">/^). It is then clear that 
a;2 = ±1 + C(e~^'^2~'^ ^Z'*) and the infimum in (7.15) is 



> hl'e 



-2S2/hi 



1 - 0{e' 



')) > I&2IV 



0{e 



-{S'+S2)/h 



)■ 



This gives Proposition 7.3 and Theorem 7.2 



7.5 An example with three minima and three saddle points 

We shall finally briefly discuss an example illustrating the generic and general cases, 
for which computations can be made by hand. We study the case of three saddle 
points and three local minima as illustrated by Figure |4j It is natural to denote 
simply by m^, fc = 1, 2, 3 the local minima and Sj, j = 1, 2, 3 the saddle points. The 



generic case corresponds to the case when Hypothesis 5.1 is satisfied. An interesting 



situation, entering in the previous general case, appears when depends smoothly 
on some parameter such that we may have 



0(mi) = 0(m2) = 0(^3), and 0(si) = 0(^2) = ^(sg). 



(7.16) 
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Figure 4: Case of 3 saddle points at the same level 



We will restrict the attention to a small neighborhood of a parameter value for which 



(7.16) holds. 



In this situation we can describe in some detail some facts about the restriction 
of the twisted Laplacian — A)^ to the space E^'^^ corresponding to the exponentially 
small eigenvalues. We know from the previous sections that this space is of dimension 
3, since this is the number of local minima. The question is to study the eigenvalues 
in more detail, and we already know that one of them is and that they are all real. 

It will be convenient to use a different system of critical components and define 
E{mk) to be the connected component of (t>^^(\ — oo, minj^^ 4'i.Sj) [) that contains m^ 
and to define Xk correspondingly as in Section |4J The eigenspace i^^"^ is generated 
by 0-forms of the form 



,(0) 



/i-"/%(/i)Xfc(x)e~('^-^(™^-))/'' + 0{e-^'^^) 



where Ck{h) ~ Ckfi + hck,i + ■■■ is a normalization constant with Ckfi > 0. Similarly we 
know that 1-forms generating the 3-dimensional eigenspace E^^^ of — A)^ associated 
with exponentially small eigenvalues can be chosen of the following form: 



where aj{x;h) ~ cijfi{x) + haj^i + ... in C°° (neigh (sj,]R")) and Vj are as in (3.12). 



Both families (e 



(0)> 



and (e\ ] 



are bases that we can assume to be orthonormal for 



the A, K scalar products 

The r 
the form 



The matrix R oi d^: i?*^*^^ — )■ E'-^^ with respect to the bases (e^ ) and (e^- ) takes 



R = h^/\pj,ke 



{sj)/h (j>{mk)/h\ 



where \pj,k\ ^ 1 for j ^ k, and = 0{l)e ch when j = k, uniformly for the parameter 
in a small neighborhood of the value where (7.16) holds. The Maxwellian e"'^/'' can 
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be written 



-<i>lh _ 7,5 \- (l + 0(e A)) ^^ (0) 



h^Y. 



-e ^ el 
^ r.u { h. ) 

k 

and is annihated by d^, so we know that the coefficent vector belongs to the kernel 
of R. It follows that 



Yl Phk- = 0{e c\), 



Cfc 
fe; k^j 



SO 



where 



P3,k=plk + 0{e-^, 



0, A; = j, 

p°fc = < djCk, k = j + 1, , rfj- X 1, Cfc X 1 
[-(ijCfc, /c = j + 2 

with the cyclic convention for the indices. Introducing 



we get 
where 



/ as -&3\ 

-Ro = -&i as , 

V ai -^2 / 

ai = o"3Ati, 02 = (Tifi2, as = ^^2/^3, 

61 = (T2/il, &2 = 0-3122, h = O-ljJ^s. 



>1^' 



We see that -Ro ( A^s ""^ =0 reflecting the fact that d^{e '^/^) = 0. This implies 
that det i?o = 0. We have 




a2 -h-A /af + b( -0162 -6103 
R*qRo = \ a2 -62 I I -bi 03 = -6201 ^2 + ^2 -«2&3 

fll -62 / \ -0361 -6302 ag + 63 

Since det(i?o-Ro) = we get that for all A G M, 

det(A - R*qRo) = \^- {al + al + al + bj + bl + bDX"^ 

+ [{ajal + bjbl + bjal) + {ajal + blbj + blaj) + {ajaj + bjbj + bjal)] A. 

(7.17) 
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To shorten the notation we let aj = Oj, f3j = b'j and 7^ = aj + (3j for j = 1,2,3. 
Then one of the eigenvalues of RqRq is and the other two are positive and given 
by 

[ '2 j -^ 

Notice that we already know that the eigenvalues are real nonnegative thanks to 



Section 3^ In particular we have 

0<D<f^^^±^^y. (7.18) 

We can localize the eigenvalues a little more: we get for example that 

n / \ ^ 7i + 72 + 73 . , ^ 

< A_ < < A+ < 7i + 72 + 73 

and also from the fact that that A+A_ = D, we get 

D , 2D 

< A_ < 



71 + 72 + 73 7i + 72 + 73 



Now we can study the case when (7.16) occurs, as illustrated by Figure 4 This 



case is included in the general situation and no longer in the generic one, since in 
particular 3 critical connected components appear when we are at the critical level 
(see Figure ^. If we assume in addition that 61 = 62 = &3 and ai = 02 = as, 
which corresponds roughly to the rotation invariant case, then we get aj = (3j =: a, 
7j = 2a for j = 1, 2, 3, and D = 9a^ = (^2i±:^±}a^ ^ gQ ^^iat the eigenvalues A+ and 
A_ of i?o are equal. 

A+ = A_ = 3a ^ ;^g-2(0(s)-0M)/h^ 

where s (respectively m) is any of the saddle points (respectively minima). 



Remark 7.4. The right-hand side inequality in (7.18) also reads 

0< (71 + 72 + 73)^ -4^) 
= (ai + a2 + as + /3i + /32 + PsY 

- 4(a2a3 + (32(3s + (32as) - ^a^ai + (3s(3i + /^aai) - 4(aia2 + /3i/32 + /3ia2)) 
= (ai - /33)' + ("2 - /3i)' + (a3 - ^2)^ 

- 2(ai - /33)(a2 - A) - 2(ai - /33)(a3 - P2) - 2(a2 - /3i)(a3 - ^2) 

Let us denote by V this last expression. From the general study we know that is 
an eigenvalue and this implies 

aia2a3 = /3i/32/33 (7.19) 

since aj = Oj, (3j = b'j, and 

det(A - RlRo) = A^ - {aj + al + al + bl + bl + bl)\^ 

+ [{44 + blbl + blal) + {alal + 6^6? + 6^a?) + {ajal + bjbl + bjal)] A 

- (010203 - b,b2b3Y. (7.20) 
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Just notice that without the assumption (7.19), V may be negative: this can happen 



def 

for example if we suppose ai — (5^ = a2 — l^i = d^ — P2 = ^ > since we get then 

V = -36^ < 0. 
Note that in this case 010203 > /3i/32/33 by direct computation so that of course 



(7.19) is not satisfied. 
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